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Abstract 

A stochastic algorithm is proposed, finding the set of intrinsic p-mean(s) associated to a prob- 
abihty measure on a compact Riemannian manifold and to p e [l,oo). It is fed sequentially 
with independent random variables (Yn)nen distributed according to i' and this is the only knowl- 
edge of v required. Furthermore the algorithm is easy to implement, because it evolves like a 
Brownian motion between the random times it jumps in direction of one of the Y^., n e N. Its 
principle is based on simulated annealing and homogenization, so that temperature and approx- 
imations schemes must be tuned up (plus a regularizing scheme if u does not admit a Holderian 
density). The analyze of the convergence is restricted to the case where the state space is a circle. 
In its principle, the proof relies on the investigation of the evolution of a time-inhomogeneous 
functional and on the corresponding spectral gap estimates due to Holley, Kusuoka and Stroock. 
But it requires new estimates on the discrepancies between the unknown instantaneous invariant 
measures and some convenient Gibbs measures. 
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1 Introduction 



The purpose of this paper is to present a stochastic algorithm finding the geometric p-means of 
probabihty measures defined on compact Riemannian manifolds, for p e [l,oo). But we are to 
analyze its convergence only in the restricted case of the circle. 

So let be given v a probability measure on M, a compact Riemannian manifold. Denote by d 
the Riemannian distance and consider for p ^ 1 the continuous mapping 

Up : M 9 X J dP{x, y) v{dy) (1) 

A global minimum of Up is called a mean of v and let Aip be their set, which is non-empty in the 
above compact setting. For p = 2 one recovers the usual notion of intrinsic mean, while for p = 1 
one gets the notion of median. 

For some applications (see for instance Pennec |12j). it may be important to find Aip or at least 
some of its elements. In practice the knowledge of v is often given by a sequence Y := {Yn)neN of 
independent random variables, identically distributed according to v. So let us present a stochastic 
algorithm using this data and enabling to find some elements of Mp. It is based on simulated 
annealing and homogenization procedures. Thus we will need respectively an inverse temperature 
evolution /3 : M+ — > M+ and an inverse speed up evolution a : M+ — > M* . Typically, they are 
respectively non-decreasing and non-increasing and we have limt->+oo ft = +Q0 and lim^^+oo on = 0, 
but we are looking for more precise conditions so that the stochastic algorithm we describe below 
finds Mp. 

Let :- {Nt)t^o be a standard Poisson process: it starts at at time and has jumps of length 1 
whose interarrival times are independent and distributed according to exponential random variables 
of parameter 1. The process is assumed to be independent from the chain Y. We define the 
speeded-up process A^(°) := [Nj:"^)t^Q via 

Vt^O, Ar(") := A^*^,, (2) 

JO as 

Consider the time-inhomogeneous Markov process X := {Xt)t^o which evolves in M in the following 
heuristic way: if T > is a jump time of A^("\ then X jumps at the same time, from Xt- to 
Xt '■= expxj,_{{p/2)l3TaTd^^'^{XT_,Yj^(a))XT_Yj^(a)). By definition the latter point is obtained by 

following during a time s := {p/2)PTO:Td^~'^{XT_,Y^{a)) the shortest geodesic leading from Xt- 

to Y^{a) at time 1 (and thus may not really correspond to an image of the exponential mapping if 

s is not small enough). The schemes a and /3 will satisfy lim(_++oo CKift = 0, so that for sufficiently 
large jump-times T, Xt will be between Xt- and Yj^(a) on the above geodesic and quite close 
to Xt-- This construction can be conducted unequivocally almost surely, because by the end of 
the description below, Xt_ will be independent of ^^(q) and the law of Xt_ will be absolutely 

continuous with respect to the Riemannian probability A. It insures that almost surely, Y^(a) is not 

in the cut-locus of Xt_ and thus the above geodesies are unique. To proceed with the construction, 
we require that between consecutive jump times (and between time and the first jump time), X 
evolves as a Brownian motion, relatively to the Riemannian structure of M (see for instance the 
book of Ikeda and Watanabe [8]) and independently of Y and A^. Very informally, the evolution 
of the algorithm X can be summarized by the equation (in the tangent bundle TM) 

Vt^O, dXt = dBt + {p/2)at(3tdP-\XT_,Y^ia))Xt-Y^^a] dN^"^ 

where {Bt)t^o would be a Brownian motion on M and where (Yj^(a))t^o should be interpreted 
as a fast auxiliary process. The law of X is then entirely determined by the initial distribution 
mo = C{Xo). More generally at any time t ^ 0, denote by mt the law of Xt. 
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We believe that the above algorithm X finds in probability at large times the set Mp of p-means, 
at least if v is sufficiently regular: 

Conjecture 1 Assume that v admits a density with respect to A and that this density is Holder 
continuous with exponent a e (0, 1]. Then there exist two constants Op > 0, depending on p ^ 1 
and a, and hp ^ 0, depending on Up and M, such that for any scheme of the form 
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where h > bp, we have for any neighborhood J\f of Mp and for any mo, 

lim mJTVl = 1 (4) 

□ 

Thus to find a element of A4p with an important probability, one should pick up the value of Xt 
for sufficiently large times t. 

The constant 6p ^ we have in view comes from the theory of simulated annealing (cf. for 
instance Holley, Kusuoka and Stroock [3]) and can be described in the following way. For any 
x,y e M, let Cx,y be the set of continuous paths C := (C(i))o^i<i going from C(0) = a; to 
C(l) = y. The elevation Up{C) of such a path C relatively to Up is defined by 

UpiC) := max UpiCit)) 
and the minimal elevation Up{x,y) between x and y is given by 



Then we consider 



Up{x,y) ■- min Up{C) 



h{Up) := max Up{x,y) — Up{x) — Up{y) + minC/p (5) 

x,yeM M 



This constant can also be seen as the largest depth of a well not encountering a fixed global 
minimum of Up. Namely, if xq e A4p, then we have 

b{Up) = maxUp{xo,y) -Up{y) (6) 

yeM 

independently of the choice of xq e A4p. 
Let us also define 

f \ _ j a , if p = 1 or p ^ 2 , , 

"^^^ |min(a,p-l) ,ifpe(l,2) 

With these notations, the main result of this paper is: 
Theorem 2 Conjecture{l\ is true for the circle M = T := M/(27rZ) with ap = a{p) and hp = h{Up). 

Let us now describe a stochastic algorithm, derived from the previous one, which should enable 
one to find the p- means of any probability measure on a compact Riemannian manifold M. 

For any x e M and k > 0, consider, on the tangent space T^M, the probability measure K^^k 
whose density with respect to the Lebesgue measure dv is proportional to {1 — k \\v\\) + (where the 
Lebesgue measure and the norm are relative to the Euclidean structure on Tx^\ Denote Kx^k 
the image by the exponential mapping at x of Kx^k- Assume next that we are given an evolution 
K : M+ 3 t 1-^ Kt e R* and consider the process Z := {Zt)t^o evolving similarly to {Xt)t^o, 
except that at the jump times T of A^("), the target Y^(a) is replaced by a point Wt sampled from 

Ky (^),KTi independently from the other variables. 
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Theorem 3 In the case M = T and p = 2, consider the schemes 

r at - {i + ty 

V t ^ 0, \ ■■= b-^ ln(l + t) 

with b > h{U2), k > and c ^ 2k + I. Then, for any neighborhood Af of M.2 CLnd for any initial 
distribution C{Zq), we get 

lim PrZteTVl = 1 

t->-+aD 

where P stands for the underlying probability. 

Still in the situation of the circle, it is possible more generally to find for any given p ^ 1 similar 
schemes (where c depends furthermore on p ^ 1) enabling to find the set of means Mp (see 
Remark [36]) . Even if v satisfies the condition of Theorem [2l it could be more advantageous 
to consider the alternative algorithm Z instead of X when the exponent a in (l3|) is too small. 
Of course we also believe that a variant of Theorem [3] should hold more generally on compact 
Riemannian manifolds M. But it seems that the geometry of M should play a role, especially 
through the behavior of the volume of small enlargements of the cut-locus of points. 

Remark 4 The schemes a, (3 and k presented above are simple examples of admissible evolutions, 
they could be replaced for instance by 

r at := Ci{r,+t)- 
Vt;>0, I ^t ■■= rMn(r2 + t) 
[ Kt = C2[r^ + tf 

where Ci, C2 > 0, ri, > 0, r2 ^ 1 and still under the conditions b > b(Up), k > and c ^ 2k + 1. 
It is possible to deduce more general conditions insuring the validity of the convergence results of 
Theorems [2] and [3] (see e.g. Proposition [27] below) . 



How to choose in practice the exponents c and k satisfying c ^ 2A; + 1 in Theorem [Sl? We note 
that the larger c, the faster a goes to zero and the faster the algorithm Z is using the data (yn)nGN- 
In compensation, k can be chosen larger, which means that v is closer to its approximation by its 
transport through the kernel K.^i^^{-) (defined before the statement of Theorem [3l for more details 
see Section 5), namely the convergence will be more precise. This is quite natural, since more 
data have been required at some fixed time. So in practice a trade-off has to be made between the 
number of i.i.d. variables distributed according to u one has at his disposal and the quality of the 
approximation of AAp. 

When 1/ is an empirical measure of the form {^ll^ 6^^)/N where the xi, I e [1, Nj, are distinct 
points of the circle, Hotz and Huckemann [5] presented an algorithm finding the means with a 
complexity of order A^^. Another algorithm was proposed by Charlier [2] valid too for means 
on the circle associated to empirical measures (the latter can next be used to approximate more 
general probability measures, but it seems this is not a very efficient method, since for each new 
point added to the empirical measure, the whole algorithm has to be started again from scratch). 
But up to our knowledge, the process of Theorem [2] is the only algorithm finding p-means for 
any p ^ 1 and for any probability measure u admitting Holderian densities, even in the restricted 
situation of the circle. 

The paper is constructed on the following plan. In next section we recall some results about 
simulated annealing which give the heuristics for the above convergence. Another alternative 
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algorithm is presented, in the same spirit as X and Z, but without jumps. In Section 3 we 
discuss about the regularity of the function Up, in terms of that of ly. It enables to see how close 
is the instantaneous invariant measure associated to the algorithm at large times t ^ to the 
Gibbs measures associated to the potential Up and to the inverse temperature (3^^- The proof 
of Theorem [2] is given in Section 4. The fifth section is devoted to the extension presented in 
Theorem [3] and the appendix deals with technicalities relative to the temporal marginal laws of 
the algorithms. 

2 Simulated annealing 

Here some results about the classical simulated annealing are reviewed. The algorithm X described 
in the introduction will then appear as a natural modification. This will also give us the opportunity 
to present another intermediate algorithm. 

Consider again M a compact Riemannian manifold and denote {•,•), V, A and A the corre- 
sponding scalar product, gradient, Laplacian operator and probability measure. Let ^7 be a given 
smooth function on M to which we associate the constant b{U) ^ defined similarly as in ([5]). We 
denote by M the set of global minima of U. 

A corresponding simulated annealing algorithm 9 := (0t)t>o associated to a measurable inverse 
temperature scheme /3 : R+ R+ is defined through the evolution equation 

V i ^ 0, dOt = dBt - y VC/(6it) dt 

It is a shorthand meaning that is a time-inhomogeneous Markov process whose generator at any 
time i ^ is L^^, where 

Lp- - l(A.-/3<VC/,V-» (8) 

Holley, Kusuoka and Stroock [3] have proven the following result 
Theorem 5 For any fixed T ^ 1, consider the inverse temperature scheme 

V t ^ 0, ft = b-^ ln(r + t) 
with b > b{U). Then for any neighborhood M of A4 and for any initial distribution C{9q), we have 

lim nOt^hf] = 1 

A crucial ingredient of the proof of this convergence are the Gibbs measures associated to the 
potential U . They are defined as the probability measures given for any /3 ^ by 

:= !2EM?M)A(<i.) (9) 

Zp 

where Zp :- ^ exp{— /3U {x)) X{dx) is the normalizing factor. 

Indeed, Holley, Kusuoka and Stroock [3] show that C{6t) and fip^. become closer and closer as t 
goes to infinity, for instance in the sense of total variation: 

^hm^||/:(eO-/^ftlltv = (10) 
Theorem [5] is then an immediate consequence of the fact that for any neighborhood J\f of M, 

lim uslMl = 1 
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The constant b(U) is critical for the behaviour (jlOp . in the sense that if we take 

V t ^ 0, pt = ln(r + t) 

with r ^ 1 and b < b(U), then there exist initial distributions C{9q) such that (jlOp is not true. 

But in general the constant b{U) is not critical for Theorem [5l the corresponding critical 
constant being, with the notations of the introduction, 

b'{U) := min maxC/(xo,y) - ^ b{U) 

xosAi yeM 

(compare with where U replaces Up and where a global minimum xq e is fixed). Note that 
it may happen that b'{U) = b(U), for instance if A4 has only one connected component. 

Another remark about Theorem [5] is that the convergence in probability of 6t for large t ^ 
toward M cannot be improved into an almost sure convergence. Denote by A the connected 
component of {x e M : U{x) ^ min^/ {7 + 6} which contains A4 (the condition b > b(U) insures 
that A4 is contained in only one connected component of the above set). Then almost surely, A is 
the limiting set of the trajectory {9t)t^Q (see [10], where the corresponding result is proven for a 
finite state space but whose proof could be extended to the setting of Theorem [5|). 

Even if we won't be interested in this paper by the analogous remarks for Conjecture [U and 
Theorem [21 they should still hold under their assumptions. 

Let us now heuristically put forward why a result such as Conjecture [T] should be true, in 
relation with Theorem [5j For simplicity of the exposition, assume that v is absolutely continuous 
with respect to A. For almost every x,y e M, there exists a unique minimal geodesic with speed 1 
leading from x to y. Denote it by {'y{x, y, t))teu, so that 7(3;, y,0) = x and 7(x, y, d{x, y)) = y. The 
process {Xt)t^Q underlying Theorem [5] is Markovian and its inhomogeneous family of generators is 
{Lat,i3t)tys0: where for any a > and /3 ^ 0, La,(s acts on functions / from C^(M) via 

VxeM, L^M]i^) - lAf{x) + - lf{j{x,y,{p/2)^adP-\x,y)))-f{x)u{dy) {11) 

2 a J 

(to simplify notations, we will try to avoid writing down explicitly the dependence on p ^ 1). The 
r.h.s. is well-defined, due to the fact that « A which implies that the cut-locus of x is negligible 
with respect to u. Furthermore Fubini's theorem enables to see that the function L^^islf] is at 
least measurable. Next remark that as a goes to 0+, we have for any / e C^{M), any x e M and 
any yeM which is not in the cut-locus of x, 

w ^ n r /(7(a;,?/, (p/2)/3ac?P-^(x,y))) - /(x) 1 « jp-v \/v7ff\ ■( n^\ 

Inn = -^pd^ x, y <V/ x , 7 x, y, > 

a^o+ a. 2 

so that for any / e C^(M) and x e Af, 

hm V;3[/](a;) = Ja/(x) + ^p f(i^'-i(x,y)<V/(x),7(x,?/,0)> K^^y) 

Recall that the potential U = UpWe are now interested in is given by ([T|) and that for almost every 
(x,2/) e M2, 

V,dP{x,y) = -pd'^'\x,y)j{x,y,0) 
(problems occur for points x in the cut-locus of y and, if p = 1, for x = y), thus 

VUp{x) = -pldP-Hx,y)j{x,y,0)iy{dy) (12) 
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It follows that or any / e C (M) and x e M, 

V /3 ^ 0, lim L^,p[f]{x) = Lp[f]{x) 

a— >0+ 

Since linif^+oo ctt = 0, it appears that at least for large times, {Xt)t^o and {9t)t^o should behave 
in a similar way. The validity of Theorem [5] for any T ^ 1 and any initial distribution C{6q) then 
suggests that Conjectured] should hold. But this rough explanation is not sufficient to understand 
the choice of the scheme (at)t>0) which will require more rigorous computations relatively to 
the corresponding homogenization property. The heuristics for Theorem [3] are similar, since the 
underlying algorithm {Zt)t^o is Markovian and its inhomogeneous family of generators (-Z^at,/?^,^^ )i>o 
satisfies 

V/eC2(M), hrn WL^^^^M - hAf]L = « 
For any a > 0, /3 ^ and k > 0, the generator La,/3,K acts on functions / e C^(M) via 
VxeM, L^,pM]i^) ■■= l^f{x) + ^j f{j{x,z,{p/2)(3adP~Hx,z)))-f{x)Ky,^{dz)u{dy) 

The previous observations suggest another possible algorithm to find the mean of a probability 
measure on M. Consider the M x M- valued inhomogeneous Markov process {Xt,Y (a) ,)i>o 

JV J +1 

where {Nf°'^)t^o was defined in ([2]) and where 



V t ^ 0, dXf = dB, 



t + {p/2)l3tdP-\Xt,Y^,^,^^)^{Xt, 0) dt (13) 

Again, up to appropriate choices of the schemes {ctt)t^o and (/3t)t>0i it can be expected that for 
any neighborhood Af of A4 and for any initial distribution C{Xq), 

lim FlXteAT] = 1 

Indeed, this can be obtained by following the line of arguments presented in [11], see [1]. 

But the main drawback of the algorithm {Xf)t^Q is that theoretically, it is asking for the 
computation of the unit vector 7(Xt; ^^(a)^-^, 0) and of the distance d{Xt,Y]^{a) ^-^), at any time 
t ^ 0. From a practical point of view, its complexity will be bad in comparison with that of the 
algorithm X := {Xt)t^o. Indeed, X is not so difficult to implement, e.g. if M is a torus. For 
simplicity, consider the case M = T, identified with (— 7r,7r], and let us construct Xf for some fixed 
t > 0. Assume we are given {Yn)neN, {cts)se[o,t]j {f^s)se[o,t] and Xq as in the introduction. We need 
furthermore two independent sequences {Tn)neN and (Ki)neNi consisting of i.i.d. random variables, 
respectively distributed according to the exponential law of parameter 1 and to the centered reduced 
Gaussian law. We begin by constructing the finite sequence iTn)nelo,N} corresponding to the jump 
times of A^^"^: let Tq :- and next by iteration, if T„, was defined, we take T„+i such that 
J^""*"^ l/ttgds = Tn+i- This is done until Tjv > t, with e N, then we change the definition of 

T]\f by imposing T/v = t. Next we consider the sequence {Xn, Xn)nelo,Nj constructed through the 
following iteration (where the variables are reduced modulo 27r): starting from Xq := Xq := Xq, if 
Xn was defined, with n e JO, A^ — 1], we consider Xn+i '■= A.« + ^/Tn+l — TnVn+i- Next we define 

p-2 



X, 



n+l 




(p/2)aT„+i/3T„+i 
(p/2)aT„+i/3r„+i 



Yn+1 ~ X„ 



(K„+i-A„+i) ,if 

p-2 ^ 

{Xn+1 - Yn+i) , otherwise 



^ vr 



Then A^r has the same law as Aj. On the contrary, the construction of Xt requires the discretiza- 
tion in time of the stochastic differential equation (jlSp and the evaluation of the corresponding 
approximation errors. 

Apart from these practical considerations, another strong motivation for this paper is the 
treatment of the jumps of the algorithms X and Z, situation which is not covered by the techniques 
of [H] (to the contrary of the jumps of the auxilliary process, which can be more easily dealt with). 
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3 Regularity issues 

From this section on, we restrict ourselves to the case of the circle. Here we investigate the 
regularity of the potential Up introduced in ([1]) and use the obtained information to evaluate how 
far are the instantaneous invariant measures of the algorithm X from the corresponding Gibbs 
measures, as well as some other preliminary bounds. 

For any x e T, we denote x' the unique point in the cut-locus of x, namely the opposite point 
x' = X + 'K. Recall that for y e T\{x'}, (7(3;, y, t))tgR denotes the unique minimal geodesic with 
speed 1 going from x to y and that 6x stands for the Dirac mass at x. 

Lemma 6 For any probability measure v on T, we have for the potential Up defined in ([Tp, in the 
distribution sense, for 2; e T, 

TT"(^\ = j pip- 1) IjdP-^iy^^) - ^pirP-^Sy^ix) u{dy) , if p > 1 
''P^''' \ 2\{5y{x)-6y>{x))u{dy) , ^f P = I 

In particular if ly admits a continuous density with respect to X, still denoted v, then we have that 



Up e C\T) and 



"(^\ = j P(P- 1) ^(2/' ^) ^(dy) - pttP ^u{x') , if p > I 

P^"^' \ {u{x)-u{x'))/7r ,^fp=l 



V X e T, U, 
Proof 

We begin by considering the case where p > 1. Furthermore, we first investigate the situation 
where v = 5y for some fixed y e T. Then Up{x) = dP{x, y) for any x e T and we have seen in ()12p 
that 

V X =H y', U'p{x) = -pdP~'^{x,y)j{x,y,0) 

By continuity of Up, this equality holds in the sense of distributions on the whole set T. To compute 
Up, consider a test function 93 e C°°(T): 

r ry+-^ ry 

If' {x)Up{x) dx = p\ ip'{x){x — y)P^^ dx — p \ ip'{x){y — x)p^^ dx 
Jt Jy •^y^T^ 

= p[^{x){x-y)P-^]y+^ -p{p-l) \'^\{x){x-y)P 

■Jy 

ry 

-pVp{x) [y - x)P''^fy_^ - p{p - 1) ip{x) [y - x)P 

= 2p'KP^^(p{y') — p{p — 1) (p{x)dP'''^{y,x) dx 

Jt 

So we get that for x e T, 

C/;(x) = p{p-l)dP-\y,x)-2p7TP-'Sy,{x) 
If p = 1, starting again from 

Vx + y', U[{x) = -j{x,y,0) (14) 

we rather get for any test function ip e C'^'(T): 

r ry+T^ ry 

(f' {x)U[{x) dx = (p'{x)dx— ip'{x)dx 

Jt Jy Jy^T^ 

= 2{ip{y') - ifiy)) 



^dx 



^dx 
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so that 



'II 



2{6y - 5yl) 



The general case of a probabihty measure v follows by integration with respect to u{dy). 

The second announced result follows from the observation that if v admits a density with respect 

to A, we can write for any x e T, 



In particular, it appears that the potential Up belongs to C'^(T), if the density u is smooth. 

Let us come back to the case of a general probability measure v on T. For any a > and /3 ^ 0, 
we are interested into the generator L^t^p defined in (jlip . Rigorously speaking, this definition is only 
valid if V is absolutely continuous. Otherwise the r.h.s. of is not well-defined for x e T belonging 
to the union of the cut-locus of the atoms of v. To get around this little inconvenience, one can 
consider for x e T, (7+(x,x + 7r,t))tgK and (7_(x,x + 'n,t))te9., the unique minimal geodesies with 
speed 1 leading from x to x + vr respectively in the anti-clockwise (namely increasing in the cover 
M of T) and clockwise direction. If y e T\{x'}, we take as before (7+(x, y, t))te]R := {^[x,y,t))t<=M. 
(7_(x, y, t))tgM. Next let be a Markov kernel from to { — ,+} and modify the definition pT]) 
by imposing that for any / e C^(T), 



where d stands for the natural derivative on T. Then the function L^^pl^f^ is at least measurable. 
But these considerations are not very relevant, since for any given measurable evolutions M+ 3 
t t— > aj e R* and M+ 9 t i— > e R4., the solutions to the martingale problems associated to the 
inhomogeneous family of generators (i>at,/3jt>o (see for instance the book of Ethier and Kurtz [3]) 
are all the same and are described in a probabilistic way as the trajectory laws of the processes X 
presented in the introduction. Indeed, this is a consequence of the absolute continuity of the heat 
kernel at any positive time (for arguments in the same spirit, see the appendix). So to simplify 
notations, we only consider the case where k{{x,y), — ) = for any x,y e T, this brought us back 
to the definition (fTTT) . where {'y{x,y,t))teM. stands for {'y+{x,y,t))teR, for any x, y e T. 

As it was mentioned for usual simulated annealing algorithms in the previous section, a tradi- 
tional approach to prove Theorem [2] would try to evaluate at any time t ^ 0, how far is C{Xt) from 
the instantaneous invariant probability fJ^at,i3ty namely that associated to L^^^^^. Unfortunately for 
any a > and /3 ^ 0, we have few informations about the invariant probability fia^p of -La,/?, even 
its existence cannot be deduced directly from the compactness of T, because the functions Lq,/3[/] 
are not necessarily continuous for / e C^(T). Indeed it will be more convenient to use the Gibbs 
distribution defined in ([9]) for /? ^ 0, where U is replaced by Up. It has the advantage to be 
explicit and easy to work with, in particular it is clear that for large /? ^ 0, /i/3 concentrates around 
A^p, the set of p- means of u. 

The remaining part of this section is mainly devoted to a quantification of what separates 
from being an invariant probability of ^o,/?, for a > and /3 ^ 0. It will become clear in the next 
section that a practical way to measure this discrepancy is through the evaluation of /i/3[(L* ^[1])^], 
where L* ^ is the dual operator of L^^p in L^(^/3) and where 1 is the constant function taking the 
value 1. Indeed, it can be seen that L* ^[1] = in ^^{iip) if and only if is invariant for L^^p- 
We will also take advantage of the computations made in this direction to provide some estimates 
on related quantities which will be helpful later on. 




V X e T, 
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Since the situation of the usual mean p = 2 is important and is simpler than the other cases, we 
first treat it in detail in the following subsection. Next we will investigate the differences appearing 
in the situation of the median. The third subsection will deal with the cases 1 < p < 2, whose 
computations are technical and not very enlightening. We will only give some indications about 
the remaining situation p e (2, oo), which is less involved. 

Some other preliminaries about the regularity of the time marginal laws of the considered 
algorithms will be treated in the appendix. They are of a more qualitative nature and will mainly 
serve to justify some computations of the next sections, in some sense they are less relevant than 
the estimates and proofs of Propositions [TOl [T4l [TSl and [20] below, which are really at the heart of 
our developments. 

3.1 Estimate of L* ^[1] in the case p = 2 

Before being more precise about the definition of -Z^* ^, we need an elementary result, where we 
will use the following notations: for y e T and 5^0, B{y, 6) stands for the open ball centered at 
y of radius 6 and for any s e M, Ty^g is the operator acting on measurable functions / defined on 
T via 

VxeT, Ty,sf{x) - f{-f{x,y,sd{x,y))) (15) 



Lemma 7 For any y e T, any s e [0, 1) and any measurable and bounded functions f, g, we have 



f gTy^sfdX = f fT, 



y-s/(l-s)9dX 



Proof 

By definition, we have 

27r gTy^sfdX = g{x)f{x + s{y-x))dx 

Jt Jy-TT 

In the r.h.s. consider the change of variables z :- sy + {1 — s)x to get that it is equal to 

1 fz-sy\ 27r r 
1 9[-. ]f{z)dz = fTy s/{i~s)fdX 

which corresponds to the announced result. 

■ 

This lemma has for consequence the next result, where T) is the subspace of L^(A) consisting 
of functions whose second derivative in the distribution sense belongs to L^(A) (or equivalently to 
L2(^^) for any /3 ^ 0). 



Lemma 8 For a > and /3 ^ such that af3 e [0, 1), the domain of the maximal extension of 

i,/3 



La,i3 on L^(/i^) is v. Furthermore the domain T>* of its dual operator L*^ a in L^(/i/3) is the space 



{/ e L (^^) : exp{—/3U2)f e 1)} and we have for any f e D* , 



1 



LUf = -exp(/3[/2)5^[exp(-/3f/2)/] 



2 



exp(/?£/2) 
a(l — a/3) 



JlB(j/,(l-a/3)7r)7'y-Q/3/(l-a/3)[exp(-^?72)/]l^((iy) - ^ 



In particular, if v admits a continuous density, then V* = T> and the above formula holds for any 
feV. 
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Proof 

With the previous definitions, we can write for any a > and /3 ^ 0, 

I 

a 



where / is the identity operator. Note furthermore that the identity operator is bounded from L^(A) 
to L^(/i^) and conversely. Thus to get the first assertion, it is sufficient to show that \Ty^aj3 t^{dy) 
is bounded from L^(A) to itself, or even only that \\Ty^a(3\\i^2(^x^^ is uniformly bounded in y e T. To 
see that this is true, consider a bounded and measurable function / and assume that a/3 e [0, 1). 
Since (Ty^a/sf)'^ = Ty^a.pf'^-, we can apply Lemma[7]with s = af3, g = 1 and / replaced by to get 
that 



^(Ty^apf) 



^dX = I /^r„ dA 



1 


1 


- a/3 




1 


1 


- a/3 




1 



3(j/,(l-^*)7r) 
3{j/,{l-s)7r) 

^ -r—n I fd\ 



dX 



Next to see that for any f,geC'^ (T) , 

j gLa,fifdfifi = ^fLlpgdfifj (16) 
where L* ^ is the operator defined in the statement of the lemma, we note that, on one hand, 

jgd^fdfip = Z^' jeM-f3U2)gd^fdX 

= J/exp(/?C/2)52[exp(-/3C/2)5]d/i/3 
and on the other hand, for any y e T, 

j gTy,ai3f dni3 = Zj^^ j exp{-l3U2)gTy^ai^f dX 



so that we can use again Lemma [71 After an additional integration with respect to v{dy), (jl6p 
follows without difficulty. To conclude, it is sufficient to see that for any / e 'L^{^p), L* ^/ e L^(//^) 
(where L* ^/ is first interpreted as a distribution) if and only if exp(— /3C/2)/ £ This is done by 
adapting the arguments given in the first part of the proof, in particular we get that 

exp(^;72) ^ r / arT^^.J^..^^ < exp(2/3osc([/2)) 



■ J lB(y,(l-a/3)7r)ry _a/3/(l-a/3)[exp(-/3C/2) • ] v{dy) 



a(l-a/3) J ^^^'^^-"^^"^^'-"'^/^^-'^^^L^ L2(A)o a2(l-a/3) 



Remark 9 By working in a similar spirit, the previous lemma, except for the expression of 

is valid for any for any a > and /5 ^ such that a/3 =|= 1. The case a/? = 1 can be different: it 

follows from 

that if V does not admit a density with respect to A which belongs to L^(A), then the domain of 
definition of L* ^, is n {/ e L^(^^) : = 0}, subspace which is not dense in L^(A) and 
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worse for our purposes, which does not contain 1. Anyway, this degenerate situation is not very 
interesting for us, because the evolutions {at)t^o and {l3t)t^o we consider satisfy atPt e (0, 1) for t 
large enough. Furthermore we will consider probability measures admitting a continuous density, 

* 

a,l/a 

a,l3 



in particular belonging to L^(A). In this case, L^i/^ and L'^i/^ admit V for natural domain, as 



in fact and L* o for any j3 ^ 0. 



For any a > and (3^0 such that a/3 e [0, 1), denote r] = a/3/(l — a/3). As seen from the previous 
lemma, a consequence of the assumption that U2 is is that for any x e T, 

L*^l(x) = Jexp(/3C/2(rE))52exp(-/3C/2(x))-- 



^ exp(/3^2(^)) 
a(l — a/3) 

z z a 



JlB(s/,(i-a/3)7r)(a;)7'j/-»?[exp(-/3f/2)](x) v^dy) 



+ ^7^-57 f exp(/3[C/2(a;) -C/2(7(x,y,-r?d(x,2/)))])z.(dy) (17) 

a(i - ap) JB(x,(l-a/3)7r) 

It appears that L* ^1 is defined and continuous if i' has a continuous density (with respect to 
A). The next result evaluates the uniform norm of this function under a little stronger regularity 
assumption. Despite it may seem quite plain, we would like to emphasize that the use of an 
estimate of L* ^1 to replace the invariant measure of by the more tractable //^ is a key to all 
the results presented in the introduction. 

Proposition 10 Assume that v admits a density with respect to A which is Holder continuous: 
there exists a e (0, 1] and A > such that 

Vx,yeT, \u{y) - i^{x)\ ^ Ad''{x,y) (18) 

Then there exists a constant C{A) > 0, only depending on A, such that for any (3^1 and 
a G (0, 1/(2/3^)), we have 

\\Llpt\\^ ^ C(^)max(a/3^a'^/3l+'^) 



Proof 

In view of the expression of L* ^1(2;) given before the statement of the proposition, we want to 
estimate for any fixed x e T, the quantity 

exp(/3[[/2(x) - [/2(7(x, y, -r]d{x, y)))]) u{dy) 

J B{x,{l~-a(5)TT) 

exp(/3[C/2(x) - U2{x - r]{y - x))]) v{dy) 

C — (1 — 0/3)77 

Lemma [6] and the continuity of the density v insure that U2 e C^(T). Furthermore, since this 
density takes the value 1 somewhere on T, we get that 

^ 2^7r" ^ 27^ A (19) 

Since U2 vanishes somewhere on T, we can deduce from this bound that ||?72llco ^ 47r^j4, but for 
A > l/(27r), it is better to use (fT2]l . which gives directly ||C/2llco ^ 

Expanding the function U2 around x, we see that for any y e (x — (1 — a/3)7r, x + (1 — a/3)7r) and 
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rj e (0, 1] (this is satisfied because the assumptions on a and /3 insure that a/5 e (0, 1/2)), we can 
find z e {x — {1 — Q;/3)7r, x + (1 — Q;/3)7r) such that 

/3[U2{x)-U2{x-r^{y-x))] = f3^U!,{x){y - x) - f3ifU^{z)^-^i^^ 

The last term can be written under the form ©^(a^/?^), where for any e > 0, Oa{^) designates a 
quantity which is bounded by K{A)€, where K{A) is a constant depending only on A (as usual O 
has a similar meaning, but with a universal constant). Note that we also have f3r]U2{x){y — x) = 
0{a/3'^). Observing that for any r, s e ffi, we can find u,v e (0, 1) such that exp(r + s) = (1 + r + 
exp(wr)/2)(l + sexp('ys)) and in conjunction with the assumption ^ 1/2, we can write that 



exp(/3[C/2(x)-C/2(x-??(y-x))]) = 1 + Pr]U^{x){y - x) + OA{a^/3*) 
Integrating this expression, we get that 

exp{/3[U2{x) - U2{j{x, y, -rj))]) v{dy) 



(20) 



I 



S(x,(l-a/3)7r) 



v[B{x,{l-aP)^)]+P^U'2{x) r 



y-xv{dy) + OA{o?P^) 



(l-o/3)7r 



Recalling that v has no atom, the first term is equal to 1 — v{B{x' ^aliir). Taking into account 
(fT2]) . we have U'2{x) = —2 ^^^^^ y — x i^{dy), so that the second term is equal to 

rx+TT rx' +aj3-K 

P'qU2{x)\ y - X u{dy) - (ir^U^ix) \ y-xu{dy) 

Jx—TT Jx'—af5-iT 

= -^m-)f + OA{a'^') 

(in the last term of the l.h.s., y — x is to be interpreted as its representative in (— vr, vr] modulo 27r). 
We can now return to (|17p and recalling the expression for U2 given in Lemma [6l we obtain that 
for any x e T, 

: ^(C/^(X))2-/3(1-K^')) 



L^pHx) 



1 

a 



^— ^— (1 - u{B{x', a/3vr) - ^([/^(x))^ + ©^(a^/j^) 
a(l — ap) \ 2 



a{l — a/3) 
+OA{a^^) 



/3 



(3' 



a 



(l-a/3)2 



(C/^(x))2 + /3 K^') 



i^{B{x',af3n) 
a/3(l - a/3) 



^1 , uiBix'aBTr) \ , ^4, 



1 - a/3 

/3 



u[x') 



a/3(l - a/3) 

u{B{x',al3TT) 



a/32 



a/3 / 1 — a/3 

1 « 9 ^ r - ^(y) + OA(a/3') 

1 - a/3 27ra/3 J2./_a/37r 



z.(x') + OA{aP^) 



The justification of the Holder continuity comes above all from the evaluation of the latter integral: 

"x'+afin 



J^x' + afiiT rx' +apTT 

v{x') - v[y) dy ^ a\ |x' 
x'—af3iT Jx'—afin 



yl" dy 



2A 



(a/37r) 



l+a 



1 + a 
^ 2A(a/37r) 
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The bound announced in the lemma fohows at once. 

■ 

To finish this subsection, let us present a related but more straightforward preliminary bound. 



Lemma 11 There exists a universal constant k > such that for any s > and (3^1 with 
j3s ^ 1/2, we have, for any y e T and f e C^(T), 

\ {Tl,[gy\{x) - gy{x)f txp{dx) ^ ks^(3^ (\{dff d^ip + \ f^ d^i^ (21) 

where T*g is the adjoint operator ofTy^g in ^.'^{fifs) and where for any fixed y sT, 

\fxeT\{y'}, gy{x) := f{x)d{x,y)j{x,y,0) 
(neglecting the cut-locus point y' of y). 

Proof 

Since the problem is clearly invariant by translation of y e T, we can work with a fixed value of y, 
the most convenient to simplify the notations being y = e M/(27rZ). Then the function g = gQ is 
given by g{x) = —xf{x) for x e (— 7r,7r). 

Due to the above assumptions, s e (0, 1/2) and we are in position to use Lemma [7| to see that for 
s e (0, 1/2) and for a.e. x e (—(1 — s)7r, (1 — s)7r), 

T:[g]{x) = ^exp(/3C/2(x))T_4exp(-/3t/2)<7](x) 

with T] := s/(l — s) and where we simplified notations by replacing Tq^ and ?o^_^ by T* and T_^. 
This observation induces us to introduce on (—(1 — s)7r, (1 — s)7r) the decomposition 

TM-g = T;[5]-^r_,[5] + ^(r_,[5]-9) + Y^5 



leading to 



where 



I 



{T:[g]{x) - g{x)f fi^idx) ^ —1-^J, + —1-^J, + -^^J, (22) 



Ji := f^' (exp(/3[f/2(x) - f/2((l + r])x)]) - l)\T^Mf l^pi^x) 

J2 := {T-r][g] - gY dfif^ 



n{l-s)n 

J3 - g dup 

J-(l-s)7r 



The simplest term to treat is J3: we just bound it above by \g^ dfijs. Recalling that g ^ vr^/^, 
we end up with a bound which goes in the direction of dHJ, due to the factor 35^/(1 - sf in (j22|) 
and the fact that /? ^ 1. 

Next we estimate the term Ji. Via the change of variable z := (1 + r])x, Lemma [7] enables to 
write it down under the form 

(1 - s) f (exp(/3[C/2((l - s)z) - U2{z)]) - ifg^z) exp(/3[[/2(z) - U2{{1 - s)z)] fip{dz) 
Jt 

= 4(1 - s) f sinh2(/3[C/2((l - s)z) - U2{z)]/2)g^iz) iip{dz) 
Jt 
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Since /5s ^ 1/2, we are assured of the bounds 

\(3[U2{{l-s)z)-U2{z)]\ ^ f3\\U^L7rs 

^ 47r2/3s 



^ 27r2 (23) 

and we deduce that 

Jl ^ 16^4 cosh2(7r2)/3252 J g2 

Again this bound is going in the direction of (121 p . 

We are thus left with the task of finding a bound on J2 and this is where the Dirichlet type 
quantity dfifs wiU be needed. Of course, its origin is to be found in the fundamental theorem 

of calculus, which enables to write for any x e (1 — s)7r). 



T-ri[9]{x) - 9{x) = -7] \ g' {{1 + 7jv)x)x dv 

Jo 



It follows that 

J2 ^ TT^T]^ '^\fs{dx) C dv{g'{il+vv)x)f (24) 

J-(l-s)7r Jo 

Recalling the definition of g, we have for any z e (— vr, vr), 

{g'{z)f ^ 2{7r\nz)f + f{z)) 

where we used again that ||C/2llco ^ /? ^ 1- Next we deduce from a computation similar 

to ()23p and from r] ^ 2s that 

^^^"^ ^ exp(4.2) 



so it appears that there exists a universal constant ki > such that 

l-il3{dx) dv {g'{{l + r]v)x)) ^ h \ dv \ X{dx) T^.riy[h]{x) 

J-(l-s)n Jo Jo J-{l-s)n 

where 

VxeT, h{x) - [if'{x)f + f\x)]fipix) 
The proof of Lemma [7] shows that for any fixed v e [0, 1], 

^ T^r,v[h]ix) X(dx) ^ — i — ( h(x)X(dx) 

J-{l-s)n 1 + vr] jj 

^ h{x)X{dx) 
Jt 

Jt Jt 

Coming back to (p^ and recalling that 77 = s/(l — s), we obtain that 

J2 ^ k2S^ ( I iff + f fdfip 
\Jt Jt 

for another universal constant /c2 > 0. This ends the proof of (|2ip . 
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3.2 Estimate of i^* ^[1] in the case p = 1 

When we are interested in finding medians, the definition (jlSp must be modified into 

VxeT, Ty^sfix) := f{7{x,y,s)) (25) 

Similarly to what we have done in Lemma [TJ we begin by computing the adjoint Ty^g of Ty^g in 
L^(A), for any fixed y e T and s e M+ small enough. 

Lemma 12 Assume that s e [0,7r/2). Then for any bounded and measurable function g, we have, 
for almost every x e T (identified with its representative in (y — vr, y + vr) /, 

Tls[g]{x) = l(y-„ + s,y-s){x)g{x -S) + t{y-s,y+s){x){g{x - s) + g{x + s)) 
+ '^iy+s,y+7T-s){x)g{x + s) 

Proof 

By definition, we have, for any bounded and measurable functions f,g, 

r riz+T 
27r gTy^sf d\ = g{x)f{x + s\gn{y-x)s)dx 

Let us first consider the integral 



ry+TT ry+TT 

g{x)f{x + sign{y — x)s)dx = g{x)f{x — s)dx 



ry+n-s 

g{x + s)f{x)dx 

Jy—s 

J-y+TT-s ny+s 
g{x + s)f{x) dx + g{x + s)f{x)dx 
y+s -^y^s 

The symmetrical computation on {y — vr, y) leads to the announced result. 

■ 

It is not difficult to adapt the proof of Lemma [H to get, with the same notations. 

Lemma 13 For a > and /? ^ such that a(3 e [0,7r), the domain of the maximal extension of 
La^fS on L^(/i^) is V. Furthermore the domain of its dual operator L"^ ^ in L^(/i^) is D* and we 
have for any f e V* , 

Ll^pf = iexp(/3C/i)52[exp(-/3C/i)/] + i fr* ,[/]Krf2/)-- 

where 

T%[/] = exp(/?C/i)T^,[exp(-/3;7i)/] 
f' 2 y^~ 

In particular, if v admits a continuous density, then V* = T> and the above formula holds for any 
feV. 

To be able to consider L* ^1, we have thus to assume that v admits a continuous density, so that 
1 e D* = D. Furthermore we obtain then that for almost every x e T, 

Ll,pm = ^{U[{x)f-^U^{x) + ^(^^T;^^[l]{x)i.{dy)-l^ 
By expanding the various terms of the r.h.s., we are to show the equivalent of Proposition llOt 
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Proposition 14 Assume that v admits a density with respect to A satisfying il8\). Then there 
exists a constant C{A) > 0, only depending on A, such that for any /3 ^ 1 and a e (0, vr/?^^), we 
have 

\\Llpl\i^ ^ C(^)max(a/3^a"/3l+") 

Proof 

From (I12p and Lemma El we deduce respectively that for all x e T, 

U[ix) = - J7(x,y,0)i/(dy) 

= — TT, x)) — X + vr)) (26) 

C/f(x) = (Kx)-Kx'))/vr (27) 
On the other hand, from Lemma [12] we get that for all s e [0,7r/2) and for almost every x e T, 

Jr*,[l](x)z.(dy) 

= i'{{x + s,x + TT — s)) exp(/3([/i(x) — Ui{x — s))) 

+i/((x - s,x + s))[exp(/3([/i(x) - Ui{x - s))) + exp(/3(C/i(x) - Ui{x + s)))] 

+i/((x — TT + s,x — s)) exp(/3(?7i (x) — Ui{x + s))) 
= v{{x,x + 7r))exp(/3(C/i(x) - Ui{x - s))) + v{{x - vr, x)) exp(/3(?7i(x) - Ui{x + s))) 

+i/((x - s,x)) exp(/3([/i(x) - C/i(x - s))) + u{{x,x + s)) exp(/3([/i (x) - Ui{x + s))) 

-v{{x' - s,x'))exp(/3(C/i(x) -Ui{x- s))) - v{{x',x' + s)) exp(/3(?7i(x) - Ui{x + s))) 

This leads us to define s = a/5/2 e (0,7r/2), so that we can decompose 



-^Ll^^pHx) = h{x,s) + l2{x,s) + h{x,s) 



with 



/i(x,s) 
h{x,s) 



1 / viix — s,x + s) , , 

— vr vix) 

IT \ s 

u{{x — s, x)) — z^((x' — s, x')) 



1 / i/((x' - s,x' + s) , 
— vr u[x ) 

IT \ S 



[exp(/3(C/i(x)-C/i(x-s)))-l] 



z^((x, X + s)) — i^((x', x' + s)) 



[exp(/3(C/i(x)-[/i(x + s)))-l] 



/3(x,s) := zy((x,x + 7r)) 



exp(/3(C/i(x) - C/i(x - s))) - 1 - spU[{x) 



+u{{x — vr, x)) 



exp(/3([/i(x) - C/i(x + s))) - 1 + s/3[/((x) 



Assumption ()18p enables to evaluate /i(x, s), because we have for any x e T and s e (0, vr/2). 



z^((x — s,x + s) , , 
vr i^(x) 



1 

2i 



i 
-f 

2» J( 



(x—s,x+s) 



u[z) — v[x) dz 



\z — x\^ dz 



{x—s,x+s) 



As'' 
1 + a 
^ As" 
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By considering the Taylor's expansion with remainder at the first order of the mapping s ^ 
exp{(3[Ui{x) — Ui{x — s)]) at s = and by taking into account ([26|) . we get for any x e T and 
(0,V(2/3)), 

\l2{x,s)\ ^ 2^^eMl3\Kl^s)l3\\U[l_^s 

^ HLo exp(/3s)/3s 
vr 

1 + ttA , , ^ ^ 

^ 2 exp(7r/2)/3s 

vr 

The term Iz{x, s) is bounded in a similar manner, rather expanding at the second order the previous 
mapping and using ([27|) to see that ^ A. 

■ 

We finish this subsection with the a variant of Lemma [11) 

Lemma 15 There exists a universal constant k > 0, such that for any s > and (3^1 with 
13s ^ 1, we have, for any f e C^{T), 

f {Ty%[9y]{x) - 9y{x)f fipidx) ^ ks' /S' ( ( {dff df,^ + ( f' df,^ 

JB{y,n-s) \J J 

where Ty ^ is the adjoint operator of Ty^g in L^(^^) and where for any fixed y e T, 

I 9y{x) ■- f{x)j{x,y,0) 

{ 9y{^) '■— '^{y-n,y-s)Lj{y+s,y+n){^)9y{^) 



V X e T\{y'}, 



Proof 

As remarked at the beginning of the proof of Lemma \TT\ it is sufficient to deal with the case 
y = 0- To simplify the notations, we remove y = from the indices, in particular we consider the 
mappings g and g defined by g{x) = -sign{x)f{x) and g{x) = t{-n-s)Lj{s,n){x)g{x). 
Taking into account that g vanishes on (— s,s), we deduce from Lemmas 1121 and 1131 that for a.e. 
X 6 ( — TT + s, TT — s), 

T:[g]{x) = exp{pU2{x))T^s[exp{-(3U2)g]{x) 
This observation leads us to consider the upper bound 

{T:[g]{x) - g{x))^ fi^idx) ^ 2Ji + 2 J2 



where 



r 

Ji := r ' (exp(/3[C/2(x) - U2{x + sign(x)s)]) - l)\T^s[g]? k^p{dx) 

J —7V + S 

riT—s 

h ■■= {T^s[g] - 9f dfi^ 

J —TT + S 

The arguments used in the proof of Lemma [TT] to deal with Ji and J2 can now be easily adapted 
(even simplified) to obtain the wanted bounds. For instance one would have noted that 

rO riT—s 

{g{x - s) - g{x)f fi/sidx) + \ {g{x + s) - g{x)f ^p{dx) 

J-TT + S Jo 



J2 
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3.3 Estimate of -L* ^[1] in the cases I < p <2 

In this situation, for any fixed y e T and s ^ 0, the definition (jlSp must be replaced by 

VxeT, Ty,sf{x) := f{j{x,y,sdP-\x,y))) 
It leads us to introduce the function z defined on (y — vr, y + vr) by 



z{x) := 



X — s{x — yY , if X e [y, y + tt) 
X + s{y — xY^"^ , if X e (y — vr, y] 



(28) 



(29) 



To study the variations of this function, by symmetry, it is sufficient to consider its restriction to 
(y, y + vr). We need the following definitions, all of them depending on y e T, s ^ and p e (1, 2): 



u+ 
u+ 
v+ 
w+ 



y + {p- l)2-ps2-p 
1 

y + s'^-p 

y- ( - 1) ^"'^ - (p - 1) '^'^ ) s 



= y + vr 



Let a{p) be the largest positive real number in (0, 1/2) such that for s e (0, cr(p)), we have u+ < 
y + iT, v+ > y — vr and w+ — y > y — v+. One checks that for s e (0,o"(p)), the function z is decreasing 
on (y,n+) and increasing on {u+,y + vr). Furthermore v+ = z{u+), w+ = z{y + vr) and n+ is the 
unique point in {u+,y + vr) such that z(n+) = y. Let us also introduce u+ the unique point in 
{u+,y + vr) such that and z{uj^) = —v+. All these definitions, as well as the symmetric notions 
with respect to {y,y), where the indices + are replaced by — , are summarized in the following 
picture (drawn by our colleague Sebastien Gadat): 



w 

+ 
V 






u u / • 
+ j/ • 


/ u 




/ V 
/ + 


w 



Figure 1: The function z 

Thus for s e (0, o"(p)), we can consider 99+ : y] — > [y, u+] and V'+ : [^^+) ^ \^+^y + ^] 
the inverses of z^ respectively restricted to [y, n+] and [n+,y + vr]. The mappings and 
are defined in a symmetrical manner on [y, f-] and These quantities were necessary to 

compute the adjoint Ty^g of Ty^^ in IL^(A), for any fixed y e T and s > small enough: 
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Lemma 16 Assume that s e [0,a{p)). Then for any bounded and measurable function g, we have, 
for almost every x eT (identified with its representative in [y — Tr,y + vr) ), 



+^{v+,y){3^)bP'-{3^)9{'4'-{x)) + ij'+{x)g{'ip+{x)) + \(p'+{x)\ g{ip+{x))] 
+l^y^^_){x)[i{j'_{x)g{^P^{x)) + V'+(x)(7(V'+(x)) + \ip'_{x)\ g{^-{x))] 

Proof 

The above formula is based on straightforward appUcations of the change of variable formula. For 
instance one can write for any bounded and measurable functions f,g defined on (y — vr, y + vr), 

9{x)f{Ty^s{.x))dx = f[z)g{ip+[z))\Lp\{z)\ dz 

■J{y,u+) •^{■"+,y) 

■ 

Since we are more interested in adjoint operators in L^(/i^), let us define for any fixed y e T, 
s e (0, (t(p)) and any bounded and measurable function / defined on (y — vr, y + vr), 

T*,[/] := exp(/3[/p)Tt,[exp(-/3f/p)/] (30) 

Then we get the equivalent of Lemmas [8] and [131 

Lemma 17 For a > and /? > such that s :- pa/3/2 e (0, cj(p)), the domain of the maximal 
extension of L^^p on is V. Furthermore the domain of its dual operator L"^ in L^(/i^) is 

D* and we have for any f e D* , 

Ll^f = iexp(/3C/p)52[exp(-/3C/p)/] + i fr*,[/]Kdy)-- 

^ (y. ] (y. 

In particular, if v admits a continuous density, then = 2? and the above formula holds for any 
feV. 

Once again, the assumption that u admits a continuous density enables us to consider L'^ ^li, which 
is given, under the conditions of the previous lemma, for almost every j; e T, by 

L;^l(rr) = ^([/;(x))2-^C/;(x) + i(Jr*^[l](x)K^y)-l) (31) 

We deduce: 

Proposition 18 Assume that v admits a density with respect to A satisfying il8\). Then there 
exists a constant C{A,p) > 0, only depending on A > and p e (1, 2), such that for any j3 ^ 1 and 
a e (0, cr(p)//3^), we have 

||L*_^1||^ ^ C{A,p)ui&-K{al3^,QP-^P^+P,a''l3^+'') 

Proof 

We first keep in mind that from (jl2p and Lemma [6l we have for all x e T, 



{x-yy-^u{dy)-\ {y-xY-^v{dy) 

Jx—n Jx 



(32) 



C/;(x) = Pip - 1) f dP-\y, x) v{dy) - pttP'Mx') (33) 
Jt 
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Taking into account (|3ip . our goal is to see how the terms l3{Up{x))'^ and —Up{x) cancel with some 
parts of the integral 



Jr*,[i](x)-iKdy) 



where s := pafi/2 e (0, (t(p)//3) a (0, cj(p)), and to bound what remains by a quantity of the form 
C {A^p){j3'^ s + l3sP^^ + s"-), for another constant C'{A,p) > 0, only depending on A > and 
(1,2). 

We decompose the domain of integration of v{dy) into six essential parts (with the convention 
that — vr ^ y — X < TT and remember that the points W-,v+, and w+ depend on y): 



Jl - 


-- {y 


e T 


y — TT < X < 




-- {y 


e T 


W- < X < v+} 


Js 


-- {y 


e T 


v+ < X < y} 


J4 - 


-- {y 


e T 


y < X < v^} 




-- {y 


e T 


V- < X < w+} 


Je 


-- {y 


e T 


< X < y + IT 



The cases of Ji and Je are the simplest to treat. For instance for Jg, we write that 
T:,,[l]{x)-liy{dy) = I 



^ f Tl,mx)-\v{dy) = -P r 

r 

Jx' 



li^idy) 



x'+wP-'-s 



S 

PTT 



X' 

p-2 



dy_ 
2tt 



ulx ) 

2 ^ ' 27rs 



r 

Jx' 



x'+nP-'^s 



v{y) - v{x') dy 



A similar computation for Ji and the use of assumption (|18p lead to the bound 



^ f Tl,[t\{x)-lv{dy)+p^^-Mx') 



vr 



^ Ap- 



(l+a)(p-l)-l 



1 + a 



^ 27rAs" 



(34) 



The most important parts correspond to J2 and J5. E.g. considering J5, which can be written 
down as the segment (x_,x+), with 



X- 

x+ := X 



X — vr + TT^ 



/ p-i 1 \ 1 

Hp - 1) 2-p _ (p _ 1) 2-p j s 2-p 

we have to evaluate the integral 

V+(^) exp(/3[C/p(x) - Up{i;+{xm - 1 u{dy) 

JX- 



p 

S 



(35) 



{y is present in the integrand through and ^'^(x)). Indeed, in view of (|32]) and ([33]) . we 

would like to compare it to 



pU'p{x) r {x - yy-^ v{dy) + p{p - 1) T {x - yf'^ u{dy) 
Jx— Jx— 



(36) 
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To do so, we will expand the terms ip'^{x) and exp(/3[C/p(x) — Up{ip+{x))]) as functions of the 
(hidden) parameter s > 0. Fix y e and recall that it amounts to x e {v-,w+). Due to ([29|) and 
to the definition of ip^, we have for such x, 



Let us begin by working heuristically, to outline why the quantities (j35|) and (j36j) should be 
close. From the above expression, we get 

By definition of ^+ , we have 

x-y = tl^+{x) - y - s{ip+{x) - y)P~'^ 

= (V^+(x)-y)(l-s(^+(x)-y)P-2) (38) 

so that X — y ^ ip^ {x) — y and 

iP'^{x) - I + s{p - l){x - yY^'^ 

On the other hand, 

exp(/3[C/p(x) - ?7p(V+(x))]) ^ 1 + /3[C/p(x) - ?7p(V+(x))] 

^ l + /3C/;(x)(x-V+(x)) 

= l-sP%[x){ij+{x)-yY-^ 

^ I - spU'j,{x){x - yf-^ 

Putting together these approximations, we end up with 

^'+(x)exp(/3[f/p(x)-f/p(^+(x))])-l - s[[p-l){x-yY-^-P%(x)[x-yY-^] 

suggesting the proximity of (f35]l and (1361) . after integration with respect to vi^dy) on (x_,x+). 

To justify and quantify these computations, we start by remarking that ip+{x) — y is bounded 
below by n+ — y, itself bounded below hy u+ —y = s'^-p . But this lower bound will not be sufficient 
in ([38|) . so let us improve it a little. By definition of u+, we have 

V- — y = u — y — s{u — yY~^ 

so that u+ — y = kpS'^-p where kp is the unique solution larger than 1 of the equation 

kp-kP-^ = (p-i)i^ -(p-i)2^ (39) 

It follows that for any y e J^, 

1 1 

1 , , , . — ^ — K ^ 



1 - s(V'+(x) - y)P-2 l-s(u+-y)P-2 



V- -y 

Kp (40) 



where the latter quantity only depends on p e (1, 2) and is given by 

K - ^ 

[p - 1) 2-P _ (p _ 1) 2-p 
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In particular, coming back to (j37p and taking into account (j38p . we get that for y e J'^, 



\ ^ I U( I ( \ ^p-2| {s{p - l)i^p+{x) - y)P 



{p-iri 



To complete this estimate, we note that in a similar way, still for y e J5, 

|(V+(x)-2/r2_(^_y)P-2| = (^_y)P-2|i_(i_,(^^(^)_y)p-2)2-p| 

^ (x-y)f-2|i_(i_,(^_^(^)_y)P-2)| 

so that in the end, 

\i>'+{x) - I - sip - 1){X - yY-^\ ^ [(p- l)2i^j2p-3)+ ;^jg2(^_y^2{p-2) ^4^-, 

We now come to the term exp(/3[J7p(x) — Up{ip+{x)y\). First we remark that 

^ pT:P-^s{ijj+{x) -yf-^ 

It follows, recalling our assumption /3s ^ <7(p), that 

|exp(^[[7p(x)-[7p(^+(x))])-l-/3[t/p(x)-t/p(V+(.x))]| ^ /3^[^p(x) -^^p(V'+(^))]% ^p(^^2^^^ 

^ 27r^/32exp(27rV(p))s2 

In addition we have, 

II [/"II 

In view of ([33]) and taking into account that J Up d\ = 0, we have 

du 



2-K 



l|c^;iL ^ 2p{p-i)M^ r 

Jo 

= 2p7r?'-^(l + vr^) 

So we get, 

^ 2^(l + ^^)(x-V+(x))' 

^ 27r(l + 7r^)s2(^^(x) -y)2(p-i) 

^ 27r^(l +7rA)s^ 
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namely 

\Up{x) - Up{^P+{x)) + sU'^{x){i,+ {x) - ^ 27r3s2 

Finally, using the inequality 

V ^ 0, Vpe (1,2), \vF'^ - v''-^\ ^ \u - v\p~^ 

it appears that 

|(^+(x)-yri-(x-yrH ^ |V+(x)-xr^ 



^ (42) 



so we can deduce that 



|exp(/3[C/p(x) - [/p(V+(:r))]) - 1 + Ps%[x)[x - yY-^\ 

From the latter bound and (|41|) . we obtain a constant K[p, j4) > depending only on p e (1, 2) 
and ^ > 0, such that 



p iP'^ix) exp(/3[C/p(x) - Up{i;+{x))]) - (1 + - l3sU;{x){x - yy-') u{dy) 

^ K{p, A) (^(3sP'^ + ph + s J''^ (x - y)2(P"2) jy{dy)j 

This leads us to upper bound 

p{x-yfP-^K{dy) ^ M^p(x-y)2(f-2)d2/ 



(43) 



with 



IX — 

1 + An r---^-'^ 
^ 1 y ^ dy 



(p- l)2-p -{p- l)2-p (44) 



An immediate computation gives, for p e (1, 2), a constant k' > such that for any s e {0,a{p)), 



tt-ttp-Is I 1 ,ifp>3/2 

r y'^^-'^dy ^ 4 J ln((l + a(p))/s) ,ifp = 3/2 (45) 

■^""p'^ [ , if P < 3/2 

Since 1 + ^ > p - 1, /3 ^ 1 and s e (0, a{p)), we can find another constant K'{p,A) > such 
that the r.h.s. of (j33]) can be replaced by K'{p, A){I3sP~^ + fi'^s). It is now easy to see that such 
an expression, up to a new change of the factor K'{p,A), bounds the difference between ([35]) and 
(j36|) . Indeed, just use that 

,,2(p-2) 



r7T — 7V^ S n7T — 7T^ 



y''^-'>dy 
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and resort to (j45p . 

There is no more difficulty in checking that the cost of replacing x_ and x+ respectively hy x — n 
and X in (IMD is also bounded by + s(p-i)/(2-p)) ^ 2K"{p, A)i3sP-\ for an 

appropriate choice of the factor K"(j), A) depending on p e (1,2) and ^ > 0. 

Symmetrical computations for J2 and remembering (j34p lead to the existence of a constant 
K"'{p, A) > 0, depending only on p e (1,2) and ^ > 0, such that for /3 ^ 1 and s e (0, a{p)/l3), we 
have 



/3(C/;(a;))2-C/;'(x) + ^Q 



T* [l](x)z.(dy)-l 



Ji u J2 u Js u Je 



^ K"'{p,A){s'' + ps- 



It remains to treat the segments J3 and J4 and again by symmetry, let us deal with J4 only: it 
is sufficient to exhibit a constant K^^'){p,A) > 0, depending on p e (1, 2) and A > 0, such that for 
^ ^ 1 and s6 (0,ct(p)//3), 



f r*,[i](x)-iKdy) 



K^^^(j),A)s^ 



(since the r.h.s. is itself bounded by K^^\p, A){a{p)) ^-p or equivalently 

T*,[l](x)-lz.(dy) 



J: 

The constant part is immediate to bound: 

J J A 



K^^Hp,A) _j_ 
^ — — -s'^-p 



P 



Idy 



1 + -kA 
27r 



r . 



/(2-p) 



Idy 



= (l + vrA)Kp 1 
27r 

For the other part, we first remark that for y e J4, we have 



(46) 



y < 



< y + HpS^-p 



y -\- g2-p < i^+ix) < y + kpS 



2-p 



2/ — (p — 1) 2-p s2-p < ip_(x) < y 

1 1 1 

y — g2-p < 'i/;_(x) < y — (p — 1) 

^ 1 

(recall that 2+ = y + kpS'^-p with /cp defined in (j39p ). It follows that we can find a constant > 0, 
depending only on p e (1, 2), such that for s e (0, o'(p)), 

max(|C/p(x) - C/p(V+(x))| ,\Up{x) - Up{i;.{x))\ , |C/j,(x) - Up{^.{x))\) ^ k''s^p 



p-i 



K';{a{p))^-ps 



In particular, we can find another constant Kp > 0, such that under the conditions that /? ^ 1 and 
Pse{0,a{p)), 

exp (/3max(|C/p(x) - ^,(V+(x))| ,\Up{x) - Up{i;^{x))\ , |C/p(x) - C/p((^_(x))|)) ^ < 
Thus, denoting one of the functions ip+, ip- or and remembering the bound ^ 1 + vrj4. 



it is sufficient to exhibit another constant > such that 



J 

Jj4 



\iIj'{x)\ dy ^ K^'^)s2-p 
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(47) 



Let us consider the case = the other functions admit a similar treatment. We begin by 
making the dependence of ip+{x) more expUcit by writing it ip-^-{x,y). From the definition of this 
quantity (see the first hne of ([38]) ) and from ([37|) . we get 

dyij+[x,y) = 



1 - s{p- l){il;+{x,y) - y)P ^ 
= -s{p-l){'il^+{x,y) -yY'^d:^ij+{x,y) 

so that the l.h.s. of (1471) can be rewritten 



— ^— I \{'P+{x,y) -yf Pdyij+{x,y)\ dy ^ ^ | {kpS^-vf P\dy^+{x,y)\ dy 
sip - 1) Jj4 ' s(p - 1) Jj^ 



k P f 

1 1 
Checking that J4 = (x — KpS ^-p , x) , the last integral is equal to 1^+ {x,x) — ip^ {x,x — KpS ^-p ) 



By 



definition of we have i/j+ix, x) = x and it appears that the quantity C '■- "P+ix, x — KpS^-p) — x 
is a positive solution to the equation 



p-i 



C = sic + KpS^-p) 

It follows that C = kpS"^ where k'p is the unique positive solution of k'p = (k'p + KpY~^ . 
Thus (j47|) is proven and we can conclude to the validity of 



To finish this subsection, here is a version of Lemma [15] for p e (1,2), which is a little weaker, since 
we need a preliminary integration with respect to viy): 

Lemma 19 Under the assumption (D3^, there exists a universal constant k{p,A) > 0, depending 
only on p e (1,2) and A > 0, such that for any s > and /3 ^ 1 with fis ^ <y{p), we have, for any 
feC\T), 



u{dy)\ {T*,,[gy\{x) - gy{x)f iip{dx) 

JT JB{y,-K-TTP-^s) 

^ k{p,A){s''^-^Up^s^)l^\^{^ffd^,p+\^f^d^^p^ (48) 
where T*,, is the adjoint operator ofTy^g in L^(yU^) and where for any fixed y e T, 



V X e T\{y'} 



9y{x) ■- f{x)dP ^{x,y)^{x,y,0) 

9y{x) ■- 1 11 Xx)gy{x) 

(y—n,y—s )u(y+s jy+ir) 



Proof 

We begin by fixing y e T and by remembering the notations of the proof of Proposition [18] (see 

Figure [TJ. Due to fact that gy vanishes on (S_,n+) = {y — s^-p ,y + s'^-p), we deduce from 
Lemma [TU] and (|30|) that for a.e. x e {y — n + 7rP~^s, y + n — ttP^^s), 

r*.a](x) = V;(x)exp(/3[[/p(x)-[/p(V'.(x))])^,(V^,(x)) 

where e e {— , +} stands for the sign of x — y with the conventions of the proof of Proposition 1181 
Thus we are led to the decomposition 

f {T;^,[gy]{x) - gy{x)f fi^idx) ^ SJi{y) + ^My) + ^My) 

JB{y,iT-TTP-^s) 
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where 



Ji{y) - \ (exp(/3[[/p(x)-i7p(Ve(:r))])-l)2(V;(x)5,(V.(x)))V/3(rfa;) 

JB(-y,7r-7rP-is) 

My) ■= {ilj'^{x)f{gy{ipe{x))-gy{x)fiip{dx) 

Jb(jj,-k--kP-^s) 

My) - f {i;U^)-l)^gl{x)f,^{dx) 



We begin by dealing with Ji{y), or rather with just half of it, by symmetry and to avoid the 
consideration of e: 

" '(exp(/3[C/p(x) - - l)2(V>^(x)5,(V+(^)))' ^dx) 

■Jy 

Let us recall that x = iIj+{x) — s(^/;+(x) — y)"^'^ and that ^p+{x) — y ^ s^^. From (|37|) we deduce 
that for X e (y, y + vr — vr^^^s), 1 ^ ^+(2;) ^ 1/(2 — p). Thus it is sufficient to bound 



r 



(exp(/3[C/p(x) - - l)^a(V+(x)))V/3(da;) 

Furthermore, for x e {y,y + -k — 7r^~^s), we have 



X-V+(x)| ^ STT^"^ (49) 



so under the assumption that s/3 e (0, 1/2), we can bound (exp(/3[C/p(x) — Up{'ip+{x))]) — 1) by a 

fyi 



term of the form kp'^s^ for a universal constant A: > 0. It remains to use g'i{x) ^ 7r^/^(x) to get 



an upper bound going in the direction of (|18|) . 

We now come to J2{y) and again only to half of it: 



ry + TT — irP 

(^;(x))2(gj^(V'+(x)) - gyix))"^ t^pidx) 

Jy 



Due to the upper bound on -0+ seen just above, it is sufficient to deal with 



ry+n—n'' 

idyii^+ix)) - gy{x))^ fi/sidx) 

Jy 



But for X e {y,y + n — n^^ s), we have ip+{x) e (y + s'^-p ,y + vr), so that gy{il)+[x)) = gy{'ilj+{x)) 
and the above expression is equal to 

"y+n—nP^^ 



ry+n—n'^ s 

{gy{i{;+{x)) - gy{x))'^ fi/sidx) 

Jy 



Coming back to the definition of gy, it appears that for x e (y,y + vr — vr^ ^s), both ip+{x) and x 
belong to the same hemicircle obtain by cutting T at y and y', so 

{gy{i;+{x)) -gy{x))'^ 
= {dP-\y,^l^+ix))f{^P+{x)) - dP-\y,x)f{x))' 

^ 2d'^P-'\y,4^+{x)){f{i;+{x))-f{x))^ + 2f\x){dP-Hy,i;+{x))-dP-Hy,x))' 
^ 27r2(P-i)(/(0+(x)) - /(x))2 + 27r2(P-i)%2{p-i)/2(^^ 
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where we have used (|42p to majorize the last term. From ()49p . we deduce that 



(/(V'+(x))-/(x))2 ^ 2s7rP-i {f{z)Ydz 



As usual, the assumption < sf3 ^ 1/2 enables to find a universal constant k > such that for 
any z e {x — sttP~^,x + sttP^^), we have Hjiix) ^ k^p[z). From the above computations it follows 
there exists another universal constant A;' > such that for any y e T, 



J2{y) ^ fc'(.^(^-i)J/2d/i^ + s2j(/')'rf/^/3 



Finally we come to Jsd/), which will need to be integrated with respect to u{dy). From (|37p . 
we first get that 

J,. f ( s{p-l)dP-^{^e{x),y) V 2, ^ . , ^ 



^ ,^ .r. S 



(2-p)2 



f d^(-p-^\Mx),y)f\x)f,^idx) 

JB(v,tt-ttP-^s) 



1 

Next, recalling that Wi^W^ ^ 1 + ttA and that d{'ip£{x),y) ^ s^-p for any x e B{y,Ti — n^^ s), it 
appears that 

I + ttA 7r2(p-i)(p- 1)2 

- s I U, (/ 



i My)u{dy) ^ ^^^L^-JPlLs' \ dy I d^^P-^\Ux).y)f\x) ^p{dx) 



l + 7r^7r2(P-i)(p-l)2 2 



27r (2-p)2 
But for any fixed z e M/(27rZ), we compute that 



fir ^ ^|d2(*'-2)(.,y)dy = 2r — i 

Jt <i{^,j/)>s^ Js2^ y2i^ 



, if p > 3/2 
^ fc^'^ ln(l/s) ,ifp = 3/2 
s 2-p , if p < 3/2 

for s e (0,1/2) and for an appropriate constant > depending only on p e (1,2). It is not 
difficult to check that as s ^ 0+, we have 

s2 , if p > 3/2 

s2(p-i) » I sHn{l/s) , Hp = 3/2 

s^s 2-p , if p < 3/2 

It follows that for any p e (1,2), we can find a constant k'{p, A) > 0, depending only on p e (1)2) 
and ^ > 0, such that 

fj3(y)K^y) ^ k'{p,A)s'(p-^') { f\x)pp{dx) 

Jt Jt 
This ends the proof of the estimate (08]). 



28 



3.4 Estimate of i^* ^[1] in the cases p> 2 

This situation is simpler than the one treated in the previous subsection and is similar to the case 
p = 2, because for y e T fixed and s ^ small enough, the mapping z defined in (j29p is injective 
when p > 2. Again for any fixed y e T and s ^ 0, the definition (llSp has to be replaced by (I28p . 
namely, 

VxeT, Ty^sf{x) := f{z{x)) (50) 

With the previous subsections in mind, the computations are quite straightforward, so we will just 
outline them. 

The first task is to determine the adjoint Ty^s of Ty^g in IL^(A). An immediate change of variable 
gives that for any s e (0, a), for any bounded and measurable function g, we have, for almost every 
X e T (identified with its representative in (y — vr, y + vr)), 

tIs[9]{x) = l(y,^(y))(x)V''(2;)5((^(x)) 

where a := T^'^^^/{p — 1) and tp : {z{y — vr), z{y + ir)) — > (y — vr, y + vr) is the inverse mapping of 
z (with the slight abuses of notation: z{y — vr) := x — tt + ttP^^s, z{y + vr) := x + tt — -rP^^s). The 
adjoint T* ,, of Ty^s in I-'^(/^/3) is still given by (I30p . As in the previous subsections, this operator is 
bounded in L^(/U^). It follows, if admits a continuous density with respect to A and at least for 
a > and /3 ^ such that s := {p/2)a/3 e [0, a), that the adjoint L* ^ of La,i3 in L^(^/3) is defined 
on T>. In particular we can consider L* ^1, which is given, for almost every x e T, by 

i;/3i(x) = ^(t/;(x))2-^^'(x) + |^(Jr*,[i](x)Kdy)-i) (51) 

From this formula we deduce: 

Proposition 20 Assume that v admits a density with respect to A satisfying il8\). Then there 
exists a constant C{A,p) > 0, only depending on A > and p > 2, such that for any (3^1 and 
as (0, cj/(p/3^)), we have 

\\Ll^pl\\^ ^ C(Ap)max(a/?^a«/3l+'') 



Proof 

The arguments are similar to those of the case J5 in the proof of Proposition [TSl but are less 
involved, because the omnipresent term 1 — s{p — l){Tp{x) — yY^^ is now easy to bound: for any 
s e [0, (t/2], we have for any y e T and x e (z(y — vr), z{y + vr)). 



^ l-{p-l)\i;{x)-y 



p~2 



s ^ 1 



In particular we have under these conditions, 



l-{p-l) lipix) -yf ^ s 



[1,2] 



Following the arguments of the previous subsection, one finds a constant K{p,A), depending only 
on p > 2 and A > 0, such that for any /3 ^ 1, s e [0,(7/(2/3)] and x e {z{y — 7r),z{y + vr)). 



il^'+{x) - 1 - (p - 1) li^+ix) - y 



exp{(3[Up{x) - [/p(^+(x))]) - 1 + ^sign(x - y)U'Jx) \x - y 



p-2 



^ K{p,A)s^ 
^ K{p,A)f]^s'^ 
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This bound enables us to approximate T*^l{x) — 1 up to a term Op^^(/3 s ) (recall that this 
designates a quantity which is bounded by an expression of the form K' {jp, A)I3'^ for a constant 
K'{p, ^) > depending on p > 2 and A > 0), by 

{[p - 1) \iIj+{x) - y\P'^ - /3sign(x - y)Up{x) \x - y\^'^]^ s 



Next we consider 



J {y eJ : xe {z{y -TT),z{y + -K))} 

= T\[x' - s7rP-\x' + sttP-^] (52) 



in order to decompose 

|^£r,%[i](x)-iKdy) = |^£r*,[i](x)-iKd2/)-|jK[^'-5vr^'-\x' + .^^-i]) (53) 

According to the previous estimate, up to a term Op^Ai/^'^s'^) the first integral is equal to 

2 ^"^^ \j ^^ "^^y^ ^) ^(^y) ~ ^^p(^) sign(2; - y)dF'^{x, y) v{dy) 

In view of (j52p . up to an additional term C?p^^(/3^s), we can replace J in the above integrals by T. 
Thus putting together (fSTI) and (f53l) with (132]) and (133]) (which are also valid here), it remains to 
estimate 



7rP~2^(x') - -v[x' - svrP-S x' + svrP-^] 



2 

and this is easily done through the assumption (llSp . 

■ 

We finish this subsection with the equivalent of Lemma [11] 

Lemma 21 For p > 2, i/iere exists a constant k{p) > 0, depending only on p > 2, such that for 
any s e (0, a), with a := 7r^~^/(p — 1), and /3 ^ 1 with (3s ^ 1, ?i;e have, for any y e T anc/ / e C^(T), 

\ iT;^s[9y]{x) - gy{x)f f,pidx) ^ k{p)s^f3'' ( \{dff d^/? + f /' 

JB(j/,7r-S7rP-l) \J J 

where Ty ^ is the adjoint operator ofTy g in L^(^^) and where for any fixed y e T, 
Vx6T\{y'}, gy{x) := f{x)dP-\x,y)^{x,y,Q) 

Proof 

We only sketch the arguments, which are just an adaptation of those of the proof of Lemma [TT] 
Again it is sufficient to deal with the case y = 0, which is removed from the notations, and 
consequently with the function g{x) = — sign(x) f{x). As seen previously in this subsection, 

we have for s e (0, cr) and x e (— vr, vr), 

Tt[9\{x) = l(_^+,^p-i^^„3^p-i)(x)exp(^[C/p(x) - Up{'ip{x))\)ilj' {x)g{'ip{x)) 
where is the inverse mapping of (— 7r,7r) 3 x ^ x — sign(x) Recall that for x e (— vr + 
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Considering the decomposition 



TMix) - g{x) 

= iexp{P[Up{x) - Up{^{x))]) - + - g{x)) + {i;'{x) - l)g{x) 

we are led, after integration with respect to l(^_T^^s'KP-^.TT-sTrP-^)ix) f^/sidx), to computations sim- 
ilar to those of Subsections 13.11 and 13.31 and indeed simpler than in the latter one, due to the 
boundedness property described in ([Mj) . 

■ 

Let us summarize the Propositions [lOl [HI [18] and [20] of the previous subsections into the 
statement: 

Proposition 22 Assume that il8\) is satisfied and for p ^ 1, consider the constant a{p) > 
defined in Then there exists two constants a{p) e (0,1/2) and C{A,p) > 0, depending only on 
the quantities inside the parentheses, such that for any a > and /3 > 1 such that af3 < cr{p), we 
have 

Despite this bound is very rough, since we have replaced an essential norm by a norm, it will 
be sufficient in the next section, when a""^^ is small, as a measure of the discrepancy between 
and the invariant measure for L^^/j. 

4 Proof of convergence 

This is the main part of the paper: we are going to prove Theorem [2] by the investigation of the 
evolution of a type functional. 

On T consider the algorithm X := {Xt)t^o described in the introduction. We require that the 
underlying probability measure z/ admits a density with respect to A which is Holder continuous: 
a e (0, 1] and ^ > are constants such that (|18p is satisfied. For the time being, the schemes 
a : M+ M* and f3 : M+ — > M+ are assumed to be respectively continuous and continuously 
differentiable. Only later on, in Proposition 1271 will we present the conditions insuring the wanted 
convergence ([4|). On the initial distribution mo, the last ingredient necessary to specify the law of 
X, no hypothesis is made. We also denote mt the law of Xt, for any t > 0. From the lemmas given 
in the appendix, we have that mt admits a density with respect to A, which is equally written 
rut- As it was mentioned in the previous section, we want to compare these temporal marginal 
laws with the corresponding instantaneous Gibbs measures, which were defined in ([9|) with respect 
to the potential Up given in ([1]). A convenient way to quantify this discrepancy is to consider the 
variance of the density of mt with respect to /i^^ under the probability measure /i^^: 



Vt>0, c;{p)It := 




Our goal here is to derive a differential inequality satisfied by this quantity, which implies its 
convergence to zero under appropriate conditions on the schemes a and /3. More precisely, our 
purpose it to obtain: 

Proposition 23 There exists two constants ci{p,A), C2{p,A) > 0, depending onp^l and A > 0, 
and a constant <;(p) e (0,1/2), depending on p ^ 1, such that for any t > with Pt ^ 1 and 
< atPt ^ ?(p), we have 

I't ^ -ci(p, A)(/3r'exp(-6(;7p)/3t) - af^^fi^ - \fi't\)It + C2{p, A){af''^ /3f + \/3't\)^/Tt 
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where b{Up) was defined in f^, a{p) in Proposition \22\ and 

a{p) := 



1 , if p = 1 or p ^ 3/2 

2{p-l) , if pe (1,3/2) 

At least formally, there is no difficulty to differentiate the quantity It with respect to the time 
t > 0. But we postpone the rigorous justification of the following computations to the end of the 
appendix, where the regularity of the temporal marginal laws is discussed in detail. Thus we get 
at any time t > 0, 

= 2 j - 1^ dtmt dX - j - dtlnifipj dfifs, - 2 j - dt\n{fi^J d^j,i3^ 
: - l) dtmt dX + 115. InML ( J - ' '^'^ + ' I 



^ 2 



dfip^ 



where we used the Cauchy-Schwarz inequality. The last term is easy to deal with: 
Lemma 24 For any t ^ 0, we have 

Proof 

Since for any t ^ we have 

VxeT, ln(^f^J = -f3tUp{x) - In (^^exp{-l3tUp{y))X{dy)^ 

it appears that 

VxeT, dtln{np^) = f^^ ^ Up{y) - Up{x) fifj^{dy) 

so that 

WdtHi^pJW^ ^ osc{Up)\f3l\ 

The bound osc(f/p) ^ vr^ is an immediate consequence of the definition ([T]) of Up and of the fact 
that the (intrinsic) diameter of T is vr. 

■ 

Denote for any t > 0, ft '■= mt/fijSf. If this function was to be C^, we would get, by the 
martingale problem satisfied by the law of X, that 
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where L^t^jSt, described in the previous section, is the instantaneous generator at time t ^ of X. 
The interest of the estimate of Proposition 1221 comes from the decomposition of the previous term 
into 

J^at,ft [ft - 1] {ft - 1) d/^ft + J^at.ft [ft - 1] rfptft 

= Jia.A [ft - 1] {ft - 1) d/i^, + J(/t - dm 

^ J [/< - 1] (ft - 1) + V^V^ft[(^*„ft[l])'] 

It follows that to prove Proposition [231 it remains to treat the first term in the above r.h.s. A first 
step is: 

Lemma 25 There exist a constant c-i{p,A) > 0, depending on p ^ 1 and A > and a constant 
o{p) e (0,1/2), such that for any a > and (3^1 such that af3'^ ^ have, for any 

/eC2(T), 

J L,,^ [/ -!](/- 1) ^ - - C3(P, ^)a"(^)/3=^) J(5/)2 d/z^ + c^ip, ^)a^(*')/33 J(/ - 1)^ dfi^ 

where a{p) is defined in Proposition \23[ 

Proof 

For any q > and /3 ^ 0, we begin by decomposing the generator L^^js into 

La,/3 = Lj3 + Ra,p (56) 

where Lp := (5^ — /3Upd)/2 was defined in ([S]) (recall that J7p is well-defined, since u has no atom) 
and where Ra,i3 is the remaining operator. An immediate integration by parts leads to 

jL[s[f-i]{f-i)dfip = -^j{^{f-l)fd^,p 

Thus our main task is to find constants cs{p,A) > and a{p) e (0, 1/2) such that for any a > 
and /3 ^ 1 with a(3^ ^ o-{p), we have, for any / e C^(T), 

Ji?„,4/-l](/-l)d/x;3 ^ csip,A)a^(-P^p'(^^{dffdii^ + ^{f-lfdt,p^ (57) 

By definition, we have for any / e C^(T) (but what follows is valid for / e C^(T)), 

VxeT, RaAf]i^) = l^f{l{x,y,ip/2)af3dP-\x,y)))-f{x),.{dy) + ^U;{x)f'ix) 

To evaluate this quantity, on one hand, recall that we have for any x e T, 

U^{x) = -p\ dP-Hx,y)j{x,y,0)u{dy) 

and on the other hand, write that for any x e T and y e T\{x}, 

f{^{x,y,{p/2)aPdP-\x,y))) - fix) = ?-a/3 Cf'{j{x,y,{p/2)af3d{x,y)u))dP-\x,y)j{x,y,0)du 

^ Jo 
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Writing s := {p/2)af3 and considering again the operators introduced in (j28p (now for any p ^ 1), 
it follows that 

jn^^p [f-i]{f-i)dfi^ 
= Y [ '^'^ J J - /'(^))(/(^) - 1)^""' 0) 

where for any fixed y e T, 

VxeT\{y}, := {f[x) - l)dP-\x,y)^{x,y,Q) (58) 

(with e.g. the convention that gy{y') '■- 0). Let us also fix the variable u e [0, 1] for a while. 

We begin by considering the case where p ^ 2. By definition of T*su (discussed in Section [3]), 
we have 

^{Ty^suifMx) - f'{x))gy{x) fifj{dx) = ^ f'{x){T*^^[gy]{x) - gy{x)) Hfsidx) (59) 

= h{y,u) + l2{y,u) 

where for any y e T, 

h{y,u) := ^ f'{x){T*^^[gy]{x) - gy{x)) fii3{dx) 

h{y,u) := - f'{x)gy{x)fip{dx) (60) 

JB(y',su'rrP-^) 

(recall from Subsections 13.11 and 13.41 that for any measurable function g, Ty ,,[g] vanishes on 
B{y' , suttP^^)). The first integral is treated through the Cauchy-Schwarz inequality, 



\h{y,u)\ ^ 



and Lemmas 1111 and 1211 at least if s/3 > is smaller than a certain constant a(p) e (0,/12). It 
follows that for a universal constant A; > 0, we have 



f \h{y,u)\ u{dy)du ^ ks^(3^ (\{df)^dfx^+ \{f-lfdfi^ C 
JTxro.i] J /Jo 



du 



<[0,1] 

2 
k 
4 



< \sp[j{^ffd^ip+\^fd^lp 

bound going in the direction of (j57p . 

Next we turn to the integral /2(y,n). We cannot deal with it uniformly over y e T but we get 
a convenient bound by integrating it with respect to ^{dy). Recalling that under the assumption 
(fT8|) the density of v with respect to A is bounded by 1 + An, it appears that 

J |l2(y,n)| Krfy) ^ l^2(y,n)|dy (61) 

^ { dy \ \f'{x)\\gy{x)\ n/sidx) 

Jt JB{y',suTTP-^) 



(1 + yl7r)7r2p-3sn f |/'||/-1| d^x^ 

JT 
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The Cauchy-Schwarz inequality and integration with respect to l[o i] (?i)(in lead again to a bound 
contributing to ([57|) . 

It is time to consider the cases where p e [1,2). We will rather decompose the l.h.s. of (j59p 
into three parts. Let us extend the notation u+ :- y + (su) ^-p from Subsection 13.31 to all p e [1, 2). 
Next we modify the definition (j58]) by introducing gy{x) := '^[y^Trfi_]u[u+ .y+w]{x)9yix) ■ Then we 
write 

j{Ty,su[f']{x) - f'{x))gy{x) fifi{dx) = Ii{y,u) + l2{y,u) + l3{y,u) 



where 



f'{^){Ty^su[9y\{x) - 9y{x)) fipidx) 

J B(y,n—sunP ^) 



h{y,u) := - f'{x)gy{x) fj.(s{dx) 

JB{y' ,suttP-^) 

Ty,su[f']ix)gy{x) fl/sidx) 



h{y,u) := 



The treatment of Ii {y, u) is similar to that of Ii {y,u), with Lemmas 1 1 5 1 and [TU] (where a preliminary 

integration with respect to i^{dy) was necessary) replacing Lemmas 1111 and 12 1[ 

Concerning /2(y,u), it is bounded in the same manner as the corresponding quantity defined in 

ioD. 

It seems that the most convenient way to deal with I^{y,u) is to first integrate it with respect to 
l[o,i](^) T-'{dy)du. Taking into account that ^ {1 + Att) and using Cauchy-Schwarz inequality, 
we get 



r X + At: f 

\ \h{y,u)\l[o^i]{u)u{y)du ^ \h{y,u)\ l[o^^{u)dydu 



1 + An 



^^ J {x) [Ty^suW] (a;))2l[o 1] {u) Hfs{dx)dydi 



1 [S_ ,u+] {x)gl (a;)l [0,1] (^) {dx)dydu 
The last factor can be rewritten under the form 



/JmMc^x) Jl^^_^^^^^^^(y)ff2(x)(iy ^ vr^-^^ J/i;3(dx)(/(x)-l)2j^ 



1 



i_ dy 



7rV2s^wJ(/-l)2o!^^ (62) 



So it remains to consider the term 

\2 



J 1 [S- ,u+] ix) {Ty^su [/'] (x) ) ^ 1 [0,1] (n) {dx)dydu (63) 
= -^^\u-,u+]{x)Ty^su[U'f]{x)iip{x)t]^Q^i^{u) dydu 

(where as a function, fip stands for the density of the measure with respect to A) . Remember that 
for any measurable function h, we have Ty^su[/i](x) := h{x + sud'^~^{x, y)'y(x, y, 0)). For x e [n_, ■U-j.], 

we have d{x, y) ^ (su) ^-p and it follows that d{x,x + sud'^ {x, y)'y{x, y,0)) ^ {su) ^-p . Taking into 
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account that ^ vr^ ^, we can then a universal constant k > such that for ^ S/S ^ (for 

an appropriate constant a{p) e (0, 1/2)) and x e T, we have Hfs{x)/fJ.i3{x + sudP~^{x,y)^{x,y,0)) ^ 
k. This leads us to consider the function h defined by 

VxeT, h{x) := {f'{x))^^ip{x) (64) 

since up to a universal constant, we have to find an upper bound of 

1 

\ T^[u-,u+]{x)Ty^su[h]{x)l[o^i]{u)dxdydu ^ \ dx \ i dy \ h{v) - 

J J-TT Jx-s'^ Jx-sdp-'^{x,y) sar {x,y) 

= \ H{v)h{v)dv 
Jt 

where for any fixed v eT, 

rjf . 1 r dxdy 

Hiv) 



S Jt2 

Let us furthermore fix x e T, 

- f«i 



^d{x,y)sCs"^ ,d{v,x)^sdP-^{x,y)^ dP ^{x,y) 



dy 2 I ( d(v, x) \ p-i 

' s - 



dP-'^{x,y) {2-p)s 
The integration of the last r.h.s. with respect to dx is bounded above by 



2 i-i^^r-'s 2 ^ 

dx = s2-p 



2-pJo 2-p 

Thus we have found a constant k{p) > depending on p e [1,2) such that (j63p is bounded above 

1 

by k{p)s^-p under our conditions on s > and (3^1. In conjunction with ()62p and definition 
(j64p . it enables to conclude to the existence of a constant k(p, A) > 0, depending on j> e [1, 2) and 
A > 0, such that 

j\h{y,u)\l[o^^{u)u{y)du ^ k{p, A)s^p^j{f - 1)2 dfi^^ J(/')2 d/x^ 

Putting together all these estimates and taking into account that /3 ^ 1, < s/? ^ a(p) and 
^2(p-i) ^ si/(2-p)^ appears that 



f 



A(y,^^) + -^2(2/,^*) + h{y,u) u{dy)du 

[0,1] 



< A)i , ifp = 1 orp^2 
^ + ,ifpe(l,2) 

< 2k'(v A)i , > ifp = 1 or 3/2 
^ 2A;(p,Aj| ^^^^2(p-i) , if (1,3/2) 



for another constant k\p^A) > 0, depending on p e [1,2) and A > 0. This finishes the proof of 

m. 



To conclude the proof of Proposition [23l we must be able to compare, for any /? ^ and any 
/ e C^{T), the energy ^pKdf)'^] and the variance Var(/, //^). This task was already done by 
HoUey, Kusuoka and Stroock [5], let us recall their result: 
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Proposition 26 Let Up be a function on a compact Riemannian manifold M of dimension 
m ^ 1. Let b(Up) ^ be the associated constant as in For any /3 ^ 0, consider the Gibbs 
measure jip given in Then there exists a constant Cm > 0, depending only on M , such that 
the following Poincare inequalities are satisfied: 

V/3^0,V/eCi(M), Var(/,/.;3) ^ Cm[1 v (/3 fij;! J]^"^-^ exp(6(^,)/3)^^[| V/p] 

We can now come back to the study of the evolution of the quantity /j = Var(/t, fip^), for t > 0. 
Indeed applying Lemma [25] and Proposition 1261 with a = at, (3 = (3t and / = ft, we get at any 
time t > such that j^t^ 1 and at/3f ^ ^(p)) 

J^ai,A [ft - 1] (/i - 

C4/3r=^exp(-6([/p)/3t) (l - 2c^{p,A)af''^(3f) Lt + c^{p, A)af''^ ti^It 



^ -(c4/3r^exp(-6([/p)A) - c^{p,A)a';'"Pt)h 
where C4 := (167r'^CT)^^ and C5(p, ^) := cz{p,A){l + 2C4). 

Taking into account Lemma [24l the computations preceding Lemma [25] and Proposition [22] one 
can find constants ci{p,A^,C2{p-,A) > and <;{p) e (0, 1/2) such that Proposition [23] is satisfied. 

This result leads immediately to conditions insuring the convergence toward of the quantity 
Lt for large times t > 0: 

Proposition 27 Let a : 1R+ M* and (3 : M+ M+ be schemes as at the beginning of this 
section and assume: 

lim /3t = +00 

>+oo 

r+ca 

{Iv ^t)-^eM-b{Up)Pt)dt = +(» 

Jo 

and that for large times t > 0, 

max{a,"(^)/3,^af(^V^M/3;|} « exp(-6(C/p)A) 
(where a{p) > and a{p) > are defined in Propositions [2E and [2^) . Then we are assured of 

lim It = 

t^+oo 

Proof 

The differential equation of Proposition [23] can be rewritten under the form 

F; ^ -VtFt + et (65) 

where for any t > 0, 

Ft 
et 



= ci(p,A)(/3r='exp(-6(C/p)/3i) -af^)/33 - |/3;|)/2 
C2{p,A){af'^(5f + m/2 



The assumptions of the above proposition imply that for t ^ large enough, /3f ^ 1 and at/Sf ^ 
?(p)) where ?(p) € (0, 1/2) is as in Proposition 1231 This insures that there exists T > such that 
([65]) is satisfied for any t^T (and also Ft < +00). We deduce that for any t ^ T, 

Ft ^ Ft exp ^- J rjs ds^ J J ^" ^'^ ^^^•^ 
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It appears that lim^^+oo -Ft = as soon as 

J -+00 
r]s ds = +00 
T 

lim et/r]t = 
The above assumptions were chosen to insure these properties. 

■ 

In particular, remarking that a{p) ^ a(p) for any p ^ 1, the schemes given in ([3]) satisfy the 
hypotheses of the previous proposition, so that under the conditions of Theorem [21 we get 

hm It = 

t—>-+co 

Let us deduce ([4]) for any neighborhood J\f of the set A4p of the global minima of Up. From 
Cauchy- Schwartz inequality we have for any t > 0, 



An equivalent definition of the total variation norm states that 

ll"^i-Mftlltv = 2inax|mt(A) -/Z;3,(A)| 
where T is the Borelian cr- algebra of T. It follows that ^ reduces to 

lim UflfTV) = 1 

/3^+oo^ 

for any neighborhood of Mp, property which is immediate from the definition Q of the Gibbs 
measures //^ for /3 ^ 0. 

Remark 28 Under mild conditions, the results of Hwang [7j enable to go further, because he 
identifies the weak limit fico of the Gibbs measures fip as /3 goes to +aD. Thus, if one knows, as 
above, that 

lim II mt — iist IL =0 

then one gets that mt also weakly converges toward fico for large times t > 0. The weight given by 



/ioo to a point x e A4p is inversely related to the value of ^Up{x) and in this respect Lemma [6] is 
useful (still assuming that v admits a continuous density). 

First note that for any x e Mp, we have Up{x) ^ 0, since x is a global minima of Up, and by 
consequence v{x') ^ 1. Next assume that we have for any x e Mp, i'{x') < 1. It follows that Mp 
is discrete and by consequence finite, since T is compact. This property was already noted by Hotz 
and Huckemann [5], among other features of intrinsic means on the circle. Then we deduce from 
Hwang [7] that 



where Z := Xja;GA4 ~ z^(x'))~^/^ is the normalizing factor. 

In this situation C{Xt) concentrates for large times t > on all the p-means of v. Thus to find 
all of them with an important probability, one should sample independently several trajectories of 
X, e.g. starting from a fixed point Xq e T. 
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□ 



Remark 29 Similarly to the approach presented for instance in [9l [TT], we could have studied 
the evolution of {Et)t>o, which are the relative entropies of the time marginal laws with respect to 
the corresponding instantaneous Gibbs measures, namely 

Vt>0, Et ■- jln(^—^ dmt 

To get a differential inequality satisfied by these functionals, the spectral gap estimate of HoUey, 
Kusuoka and Stroock [1] recalled in Proposition must be replaced by the corresponding loga- 
rithmic Sobolev constant estimate, which is proven in the same article [1]. 

□ 



5 Extension to all probability measures u 

Our main task here is to adapt the computations of the two previous sections in order to prove 
Theorem [3l As in the statement of this result, it is better for simplicity of the exposition to restrict 
ourselves to the important and illustrative case p = 2, the general situation will be alluded to in 
the last remark of this section. 

We begin by remarking that the algorithm Z described in the introduction evolves similarly to 
the process X, if we allow the probability measure to depend on time. More precisely, for any 
K > 0, consider the probability measure u^, given by 

VzgM, u^{dz) := ^ u{dy)Ky^^{dz) (67) 

where the kernel on M, {y,dz) i— > Ky^^^{dz) was defined before the statement of Theorem [31 For 
a > 0, /3 ^ and k > 0, let us denote by Lq,,/?.^ the generator defined in pT]) . where v is replaced 
by u^. Then the law of Z is solution of the time-inhomogeneous martingale problem associated to 
the family of generators {Lat,i3t,Kt)t^o- This observation leads us to introduce the potentials 

Vk>0, VxeM, C/2,«(x) ^ d^{x,y) i^^{dy) 

as well as the associated Gibbs measures: 

V/3^0, Vk>0, fi^,^{dx) := Z"]^ exp(-/3C/2,«(x)) A(dx) 

where Z^^^ is the renormalization constant. 

Denote by nit the law of Zt for any t ^ 0. The proof of Theorem [3] is then similar to that of 
Theorem [2] and relies on the investigation of the evolution of 



V t > 0, 




which play the role of the quantities defined in (|55|) . 

While the above program was presented for a general compact Riemannian manifold M, we 
again restrict ourselves to the situation M = T. 

We first need some estimates on the probability measures z/^, for k > 0: 
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Lemma 30 For any k > 0, z^^ admits a density with respect to A, still denoted z^^- Furthermore 
we have, for any k > l/n, 

where dv^^ stands for the weak derivative (so that the last norm is the essential supremum 
norm with respect to X). 

Proof 

When M = T, for any k > 0, the kernel K.^f^(-) corresponds to the rohing around T of the kernel 
defined on R by (y, dz) ^ k{1 — k\z — y\)+ dz. In particular for any y e T, Ky^f^{-) is absolutely 
continuous with respect to A and (I67p shows that the same is true for u^. If furthermore k > l/vr, 
from this definition we can write for any z e T, 



(rz+l/K 
J ^ {'^ - Kd{y,z))+ u{dy) ] dz 

namely, almost everywhere with respect to X{dz), 

rz + l/ti 

Uf,{z) = IttkI {1 - Kd{y,z))+iy{dy) 

Jz-1/k 
rz + l/K 

^ Ittk ^{dy) 

Jz-I/k 

^ 2ttk 

Next for almost every x, y e T, we have 

\v^{x) - Vf^[y)\ ^ 27rK \{l - Kd{x, z))+ - {I - K,d{y, z))+\ u{dz) 

Jt 

^ 27rK |1 — Kd{x, z) — 1 + Kd{y, z)\ v{dz) 
Jt 

^ 2t:k^ \d{x , z) — d{y , z)\ v{dz) 
Jt 

^ 2'KK^d{x, y) 

This proves the second bound. 



An immediate consequence of the last bound is that for any x e T, the map (l/vr, +oo) 9 k > C/2,k(^) 
is weakly differentiable and for almost every k > l/vr, |(?kJ72,k(2;)| ^ 27r^K^. But one can do better: 

Lemma 31 For any x eT and any k > l/vr, we have 

|5«C/2,^(x)| ^ 

Proof 

It is better to come back to the definition of to get, for x e T and k > l/vr (where stands 
for weak derivative): 

dKU2,K{x) = d,^ (^27TK j X{dy) d'^{x,y) j {1 - Kd{y,z))+ iy{dz) 

= 2-K X{dy) d^ {x , y) u{dz) {1 — Kd{y, z))^ —2-kk X{dy)(f{x,y) u{dz)d{y,z) 
J Jt J -Jy-i/n 
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The first term of the r.h.s. is equal to U2^k{x)/k, and is bounded by ||C/2,k|Ioo ^ '^'^/i^- In absolute 
value, the second term can be written under the form 



X{dy)(f{x,y)d{y,z) ^ 2TT^K\u{dz)\ X{dy) \y - z\ 

z—1/k J Jz~1/k 

27r3 



The improvement of the estimate of the previous lemma with respect to the one given before its 
statement is important for us, since it enables to obtain that if [Ptjt^o and {nt)t^o are schemes, 
then we have 

Vt^O, ||5tln(^^,,,J||^ ^ 7r2|/5;| +37r3A|(ln(Att))'| (69) 

This bound replaces that of Lemma [241 in the present context. Note that for the schemes we 
have in mind and up to mild logarithmic corrections, we recover a bound of order 1/(1 + t) for 
||(5( ln(/i^j Kt)llco' which is compatible with our purposes. 

In the same spirit, even if this cannot be deduced directly from Lemma [311 we have 

Lemma 32 As k goes to infinity, U2.K converges uniformly toward U2- In particular, ifh{-) is the 
functional defined in then we have 

lim h{U2,n) = KU2) 



Proof 

Since ||5[/2,k|Ioc ^ k > 0, it appears that {U2,k)k>q is an equicontinuous family of 

mappings. It is besides clear that i^^ weakly converges toward v as k goes to infinity, so that 
U2.K,{x) converges toward U2{x) for any fixed x e T. Compactness of T and Arzela-Ascoli theorem 
then enable to conclude to the uniform of f/2,K toward U2 as k goes to infinity. The second assertion 
of the lemma is an immediate consequence of this convergence. 

■ 

Consider for the evolution of the inverse temperature the scheme 

V t 0, A := ln(l + t) 

where h > b(U2) and denote p :- {1 + b{U2)/h)/2 < 1. Assume that the scheme {nt)t^Q is such that 
lim(_>+oc /^i = +00. Then from the above lemma and Proposition [26] (recall that ||5C/2,k|Ioc ^ 
for any k > 0), there exists a time T > such that for any t ^ T, 

V/£Ci(T), ^^^^Var(/,M^„,J ^ P^.Ai^ff] (70) 

Like (|69p , this crucial estimate for the investigation of the evolution of the quantities (j68p still does 
not explain the requirement that k e (0, 1/2) in Theorem [3l Its justification comes from the next 
result, which replaces Proposition [10] in the present situation. 

Proposition 33 For a > 0, /3 ^ and k > 0, let L'^ ^ ^ he the adjoint operator of L^^p^n in 

L^(//^.fj)- There exists a constant Ci > such that for any /3 ^ 1, k ^ 1 and a e (0,(2/5)^^ a 
(/33(/3 + k))-i/2)^ u;e have 
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Proof 



It is sufficient to replace U2 by f72,K in the proofs of Section [Sj in particular note that (I17p still 
holds. From Lemma [6] and the first part of Lemma [30l it appears that ()19p has to be replaced by 



V K ^ 1, 



\\UnJ ^ 4ttk 



Instead of ()20p . we deduce that for any x,y e T and a, f3 and k as in the statement of the 
proposition, 



exp ( (3 



a/3 



1 - a/3 



{y - x) 



1 



1 - a/3 



U^,Ax){v-x) + 0{a^/3\p + K)) 



Keeping following the computations of the same proof, we end up with 

Llp^^t{x) = ^ ^ / ^ r ^"""^ y.{x') - v^{y) dy + 0(a/33(/3 + n)) 
1- a/3 lTTa/3 J^./.^^^ 

To estimate the last integral, we resort to the second part of Lemma [30) we get 



rx' +al3iT 
Jx'—aBn 



[y) dy 



\x' — y\ dy 

x'—a/UTT 

= iTTK^ia^Trf 



This leads to the announced bound. 



Similar arguments transform Lemma [251 into: 

Lemma 34 There exists a constant C2 > 0, such that for any a > 0, /3 ^ 1 and k ^ 1 with 
a/32 ^ ^^2, we have, for any f e C^iT), 

jL,,;3,4/-l](/-l)d/i^ ^ -(^^-C2a/3\/3 + K)^ ^{dffdfip 



+C2a0- 



\f3 + K)\^{f-lfd^l^ 



Proof 

The modifications with respect to the proof of Lemma [251 are very limited: one just needs to take 
into account the bounds lit/' . 11 ^ 2tt and lli^K-ll-r, ^ 27rK for k ^ 1. Indeed, there are two main 
changes: 

• in (j56p . where the remaining operator has to be defined by 

• in (j6ip . the factor 1 + An must be replaced by 2-11 k, by virtue of the first estimate of Lemma [30l 
It leads to the supplementary term a^'^K in the bound of the above lemma. 

■ 

All the ingredients are collected together to get a differential inequality satisfied by {If)^Q. 
More precisely, under the requirement that (j70p is true for t ^ T > 0, as well as /3t ^ 1, ^ 1 
and at/3^y/Kt ^ 1/2, we get that there exists a constant C3 > such that 



Vt^T, Z; ^ -7]tlt + et^/It 
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where for any t ^ T, 

m ■■= ^^^^-Csiat/3?{(3t + Kt) + \f3[\+/3t\{ln{Kt)y\) 

et := CsiatP!{f3^ + K^t) + Wt\ + f3t\{HKt)y\) 

Under the assumptions of Theorem [3] (ah'eady partiaUy used to insure the vahdity of (j70p for some 
p e (0, 1)), it appears that as t goes to infinity, 



1 



€t = O 
and this is sufficient to insure that 



hm It 

t— >+00 



l + t)p 
1 



l + t 



The proof of Theorem [3] finishes by the arguments given at the end of Section [H 

Remark 35 As it was mentioned at the end of the introduction, if one does not want to waste 
rapidly the sample (l'n)neN (especially if it is not infinite ...), one should take the exponent c the 
smallest possible. From our assumptions, we necessarily have c > 1. But the limit case c = 1 can 
be attained: the above proof shows that the convergence of Theorem [3] is also valid for the schemes 



6-1 ln(l + t) 
ln(2 + t) 



The drawback is that v is not rapidly approached by i/^t as t goes to infinity and this may slow 
down the convergence of the algorithm toward M. Indeed, from the previous computations, it 
appears that the law of is rather close to the set of global minima of C/2,Kt ■ 



Remark 36 The cases p = 1 and p ^ 2 can be treated in the same manner, but for p e (1, 2), one 
must follow the dependence on A of the constants in the proof of Lemma [191 In the end it only 
leads to supplementary factors of k, so that Theorem [3] is satisfied with a sufficiently large constant 
c, depending on p ^ 1 and on the exponent k entering in the definition of the scheme (Kt)t>o- But 
before going further in the direction of this generalization, it would be more rewarding to first 
check if the dependence on p of Op in Theorem [2] is just technical or really necessary. 

□ 



A Regularity of temporal marginal laws 

Our goal is to see that at positive times, the marginal laws of the considered algorithms are 
absolutely continuous and that if furthermore « A, then the corresponding densities belong to 
C^{T). We will also check that this is sufficient to justify the computations made in Section 4. 

Let X be the process described in the introduction, for simplicity on T, but the following 
arguments could be extended to general connected and compact Riemannian manifolds. We are 
going to use the probabilistic construction of X to obtain regularity results on mt, which as usual 
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stands for the law of Xf, for any t ^ 0. So for fixed t > 0, let be the largest jump time of 
A^(") in the interval [0, t], with the convention that Tt = if there is no jump time in this interval. 
Denote by the law of (T^, XxJ on [0, t] x T. Furthermore, let Ps{x, dy) be the law at time s ^ 
of the Brownian motion on T, starting at x e T. From the construction given in the introduction, 
we have for any t > 0, 

mt{dx) = S,t{ds,dz) Pt-s{z,dx) (71) 

J[0,t]xT 

An immediate consequence is: 

Lemma 37 Let t > be fixed. About the measurable evolutions a : R+ R* and (3 : R+ — > 
R+, only assume that inf^gp.t] > 0- Then, whatever the probability measure v entering in the 
definition of X, we have that nit is absolutely continuous. 

Proof 

By the hypothesis on a, is the unique atom of ^(•, T), the distribution of (its mass is .^t({0}, T) = 
exp(— 5p l/cKs ds)) and C(-, T) admits a bounded density on (0, t]. Since furthermore for any s > 
and z e T, Ps{z, •) is absolutely continuous, the same is true for mt due to ()7ip . 

■ 

To go further, we need to strengthen the assumption on v. 

Lemma 38 In addition to the hypotheses of the previous lemma, assume that v admits a bounded 
density and that i^ise[o,t] Ps > 0- Then for any t > 0, the density of mt belongs to C^{T). 

Proof 

We begin by recalling a few bounds on the heat kernels Ps{x,dy), for s > and x e T. We have 
already mentioned they admit a density, namely they can be written under the form ps{x,y) dy. 
Since the Brownian motion on T is just the rolling up of the usual Brownian motion on R, we have 
for any x e T, 

V,e(x-.,x + vr], ,.(x,,) ^ ^. ''''''^-iy - ^^^n)Vi2s)) ^^^^ 

From a general bound due to Hsu [6], we deduce that there exists a constant Co > such that for 
any s > and y e (x — tt, x + vr], we have 

\dyps{x,y)\ ^ Co(^^ + ;i^)p.(x,y) 

To get an upper bound on ps{x,y) = Ps{0,y — x), consider separately in (f72]l the sums of n e 
and n e Z_o-\{0}, where a e { — , +} is the sign of y — x. It appears that for s e (0, t], 

I N ^ o V exp(-(y -x + 27rn)7(2g)) 
Ps[x,y) ^ 2 



V27rs 

^2 , 



exp(-d2(x,y)/(2s)) 



^ Ci{t) 



/27rs 



where Ci{t) := X!nGZ+ exp(— 2(7rn)^/i). Taking into account ([7T]) and Lemma[371 if we were allowed 
to differentiate under the sign integral, we would get for any x e T, 

dxmt{x) = ^t{ds,dz) dxPt-s{z,x) (73) 

J[0,i]xT 
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(where the l.h.s. stands for the density of rrit with respect to 27rA). Unfortunately the usual 
conditions don't apply here, so it is better to consider the approximation of the density mt by m^,*, 
where for e e (0, t), 

V X e T, mtAx) '■= ^t{ds,dz) pt-s{z,x) 

J[0,t-e]xT 

There is no difficulty in differentiating this expression under the sign sum and in the end it appears 
to be smooth in x. So to get the announced result, it is sufficient to see that dxm^^tix) converges 
to the r.h.s. of (173^ . uniformly in x e T as e goes to 0+. Let us prove the stronger convergence 



lim sup it{ds,dz) \dxPt-s{z,x)\ = 

^-^0+ xgT J\t-e,t\xT 



z l\ exp(-z7(2s)) 



l[t-e,t]xJ 

The assumptions that inf^g^.t] ctsPs > and that z/ admits a bounded density imply that the latter 
is equally true for ^((s, •), the regular conditional law of Xt^ knowing that Tt = s, for any s > 0. 
We can even find C2{t) > such that (,tis,dz) ^ C2{t)dz, uniformly over s e {0,t] (but a priori 
C2{t) may depend on t > through inf^gp,*] Oisl^s)- In tlie proof of Lemma [371 we have already 
noticed that there exists C3{t) > such that (,t{ds,T) ^ C3{t)ds, for s =1= 0. It follows that for 
ee (0,t), 

J it{ds,dz) \dj,pt-s{z,x)\ 

= 2CoCi{t)C2{t)C^{t) [ dz [ ds . > , 

Jo Jo \s'^' sj V27r 

This bound no longer depends on x and to compute the latter integral, consider the change of 
variable u = z^/s, z being fixed: 

{ dz{ ds ( + exp(-z2/(2s)) = { dz { du ( + u] exp(-n/2) 

Jo Jo \S-^/^ S) Jo J^2/, VV^ / 

We conclude by remarking that by the dominated convergence theorem, the latter term goes to 
zero with e. 



Remark 39 More generally, but still under the assumption that v admits a bounded density, 
the density rat is at some time t > 0, if we can find e e (0, i) such that inf^g[^_^ > and 
infgg|-f_£^^] > 0. This comes from the above proof or can be deduced directly from Lemma [38l 
and the Markov property of X. 

□ 

The same arguments cannot be used to prove that for t > 0, the density of mt belongs to C^(T). 
A priori, this is annoying, since in Section HI to study the evolution of the quantity It defined in 
()55p . we had to differentiate it with respect to t > and the computations were justified only if 
the densities mt were C^. The classical way go around this apparent difficulty is to use a mollifier. 

Let p be a smooth nonnegative function on M whose support is included in [—1, 1] and satisfying 
l^p{y) dy = 1. For any 6 e (0, 1), define 
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(where functions on T are naturally identified with 27r-periodic functions on M) . These functions are 
smooth and what is even more important for Section 4, the mapping M* x T 3 [t, x) ^ d'^rnf\x) 
is continuous. Furthermore, the mf^ are densities of probability measures on T. More precisely, 
for any t ^ 0, mf^ is the density of C{Xt) when C{Xq) = rri^\ as a consequence of the linearity 
of the underlying evolution equation (i.e. V f ^ 0, dtrrit = mtLat^pt, in the sense of distributions). 
Thus the computations of Section 4 are justified if we replace there {mt)t>o by {ml^^)t>o, for any 
fixed 6 e (0,1). In particular the inequality (j66]) is satisfied for (m^''^)f>o instead of {mt)t>o- It 
remains to let 5 go to 0+ to see that the same bound is true for the flow {mt)t>o- This proves 
Theorem [2] for general initial distributions mo, for instance Dirac masses. In fact, one could pass 
to the limit 5 — > 0+ before (|66p . for instance already in Proposition 1231 to see that it is also valid. 
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